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ABSTRACT 


The  object  of  this  study  is  to  present  new  results  concerning  the 
integral -equation-Green-function  method  of  solution  of  the  basic  linearized 
problem  of  potential  flow  of  regular  water  waves  about  a  body.  Specifically, 
this  basic  potential-flow  problem  consists  in  determining  the  velocity  potential 
given  by  the  solution  of  the  Laplace  equation  subject  to  the  usual  linearized 
boundary  condition  on  the  mean  free  surface  and  to  a  Neumann  condition  on 
the  mean  body  surface.  The  main  results  of  the  study  are:  (i)  a  new  integral 
equation  for  determining  the  velocity  potential,  (ii)  a  relatively  simple 
explicit  approximate  solution  of  this  integral  equation,  and  (iii)  asymptotic 
expansions,  ascending  series,  and  one-dimensional  Taylor  series  expansions 
for  efficient  numerical  evaluation  of  the  Green  function. 
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Introduction 

The  object  of  this  study  is  to  present  new  results  concerning  the 
integral-equation-Green-function  method  of  solution  of  the  basic  linearized 
problem  of  potential  flow  of  regular  (time-harmonic)  water  waves  about  a 
body.  Particular  problems  of  practical  interest  encompassed  in  this  basic 
potential-flow  problem  are  the  usual  problems  of  wave  diffraction,  in  which 
the  body  is  held  fixed,  and  of  wave  radiation,  in  which  the  body  is  forced 
to  oscillate  about  a  mean  position  in  otherwise  calm  water;  the  problem  of 
linearized  motion  of  a  freely-floating  body  in  regular  waves  can  be  decomposed 
into  such  a  wave-diffraction  and  six  fundamental  wave-radiation  problems 
(corresponding  to  the  six  possible  degrees  of  freedom  of  motion  of  an  unre¬ 
strained  rigid  body),  as  is  well  known.  However,  this  study  is  not  concerned 
with  these  important  particular  problems  per  se,  but  with  the  development  of 
a  general  method  of  solution  of  the  basic  potential-flow  problem  common  to 
these  practical  problems.  Specifically,  this  basic  potential-flow  problem 
consists  in  determining  the  velocity  potential  (<(>)  given  by  the  solution  of 
the  Laplace  equation  subject  to  the  usual  linearized  boundary  condition  on 
the  mean  free  surface  and  to  a  Neumann  condition  (specified  normal  derivative 
<f  )  on  the  mean  body  surface. 

The  main  results  obtained  in  this  study  are  a  new  integral  equation 
for  determining  the  velocity  potential,  and  asymptotic  expansions,  ascending 
series,  and  one-dimensional  Taylor  series  expansions  for  fast  and  easy  numerical 
evaluation  of  the  Green  function  (i.e.  the  velocity  potential  of  the  flow 
caused  by  a  source  of  pulsating  strength  at  a  fixed  position  in  otherwise  calm 
water);  furthermore,  a  relatively  simple  explicit  approximation  for  the  velocity 
potential  is  given,  although  the  practical  usefulness  of  this  explicit 
approximate  solution  remains  to  be  tested  numerically.  The  plan  and  content 
of  the  study  are  briefly  described  below. 

The  basic  potential-flow  problem  is  formulated  in  section  1.  In  section 
2,  the  Green  function  associated  with  the  linearized  free-surface  boundary 
condition  is  defined,  and  the  classical  expression  —  in  terms  of  a  double 
integral  —  for  this  fundamental  function  is  briefly  derived.  The  main  result 
of  section  2  is  stated  by  equations  (2.13)  and  (2.14).  Equations  (2.13)  are 
well  known;  these  equations,  however,  are  valid  only  for  x,  strictly  negative, 
and  equations  (2.14)  are  proper  in  the  limiting  case  r,  =  0. 
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This  distinction  is  important  in  section  3,  where  fundamental  integral 

identities  verified  by  the  velocity  potential  are  obtained  by  using  equations 

(2.13)  and  (2.14)  in  the  classical  Green  identity  (3.1).  Identities  (3.6a,b,c) 

are  well  known  [although  the  particular  case  of  these  equations  corresponding 
2  2 

to  V  <|>  =  0  =  (j>z-f(l+ie)  <{>  is  usually  given  in  the  literature,  and  the  traditional 
derivation  of  these  identities  from  the  Green  identity  (3.1)  is  somewhat  different 
from  that  shown  in  this  study].  The  main  result  of  section  3,  however,  is  the 
integral  identity  (3.9).  This  new  integral  identity  is  valid  outside,  inside, 
and  on  the  body,  and  thus  is  equivalent  to  the  set  of  the  three  classical  identities 
(3.6a,b,c).  An  interesting  alternative  form  of  the  integral  identity  (3.9)  is 
given  by  equation  (3.10),  which  involves  a  "modified  Green  function".  Although 
the  present  study  is  primarily  concerned  with  the  "exterior  potential-flow 
problem",  that  is  the  problem  of  flow  about  a  body,  the  integral  identities 
corresponding  to  identities  (3.6),  (3.9)  and  (3.10)  are  also  given  for  the 
"interior  potential-flow  problem".  It  is  interesting  that  whereas  the  classical 
identities  (3.6a,b,c)  and  (3.11a,b,c)  for  the  exterior  and  interior  problems, 
respectively,  are  identical  (provided  only  one  properly  interchanges  the  terms 
"exterior"  and  "interior"),  identities  (3.9)  and  (3.12),  or  the  equivalent 
alternative  identities  (3.10)  and  (3.13),  are  not  identical.  In  particular 
identities  (3.10)  and  (3.13),  although  identical  in  form,  involve  the  different 
"modified  Green  functions"  &  and  G'*'  defined  by  equations  (3.10a)  and  (3.13a) 
for  the  "exterior"  and  "interior"  problems,  respectively.  The  zero-and  infinite- 
frequency  limits  finally  are  briefly  considered. 

Section  4  is  mainly  concerned  with  the  usual  linearized  problem  of 
flow  of  regular  waves  about  a  body.  In  this  important  special  case,  the 
general  integral  identity  (3.9)  becomes  equation  (4.2),  which  provides  an 
integral  equation  for  determining  the  potential  <j>  on  the  body  surface.  This 
new  integral  equation  is  one  of  the  major  results  of  the  present  study,  and  it 
is  discussed  in  soi.  detail.  In  particular,  a  notable  feature  distinguishing 
the  integral  equat1  l  (4.2)  from  the  classical  integral  equation  (3.1)  and 
analogous  integral  equations  is  that  every  term  in  equation  (4.2)  is  continuous 
at  the  body  surface.  Another  interesting  feature  of  the  integral  equation 
(4.2)  is  the  presence  of  the  "waterplane  integral"  defined  by  formula  (3.9a) 
in  the  case  of  a  free-surface  piercing  body.  (This  waterplane  integral  is 
not  present  in  the  case  of  a  fully-submerged  body).  A  simple  explanation  for 
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the  occurence  of  the  "waterplane  integral"  is  provided  by  considering  a 
free-surface  piercing  body  as  the  "zero-submergence  limit"  of  a  slightly- 
submerged  body  consisting  of  the  mean  wetted  surface  of  the  body  closed  by  a 
horizontal  "lid".  A  relatively  simple  explicit  approximation  for  the  velocity 
potential  <f>,  defined  by  formulas  (4.8)  and  (4.3),  is  finally  given;  this 
explicit  approximate  solution  is  one  of  the  main  results  of  the  present  study. 

A  major  difficulty  in  the  numerical  solution  of  the  integral  equation 
(4.2)  or  in  the  evaluation  of  the  approximation  (4.8)  stems  from  the  fact 
that  calculation  of  the  basic  double-integral  representation  of  the  Green 
function  given  by  expression  (2.10b),  or  by  the  alternative  expression  (2.10a), 
requires  excessive  computing  time.  The  object  of  the  remaining  five  sections 
of  the  study  thus  is  to  obtain  simplified  expressions  for  the  Green  function 
that  can  easily  and  efficiently  be  evaluated  numerically.  In  a  first  step, 
three  alternative  and  complementary  single-integral  representations  of  the 
Green  function,  that  is  expressions  involving  single  (one-fold)  integrals 
only,  are  obtained  from  the  two  above-mentioned  double-integral  representations. 
These  three  single-integral  representations,  given  in  section  5,  are:  (i)  the 
"Haskind  integral  representation"  (5.8c),  which  is  essentially  identical  to 
an  expression  obtained  (in  a  different  manner)  by  Haskind  [1] ,  (ii)  the  "near¬ 
field  integral  representation"  (5.11),  which  was  also  obtained  (independently 
and  in  a  different  manner)  by  Martin  [2],  and  (iii)  the  far-iield  integral 
representation"  (5.21),  which  does  not  appear  to  have  been  given  previously. 

The  "near-field  representation"  (5.11)  and  the  "far-field  representation" 

(5.21)  are  analogous  to  the  single- integral  representations  of  the  Green 
function  of  ship-wave-resistance  theory  given  in  [3]. 

In  a  second  step,  asymptotic  expansions,  ascending  series,  and  one¬ 
dimensional  Taylor  series  expansions  are  obtained  from  the  three  above-mentioned 
alternative  single-integral  representations  of  the  Green  function.  Specifically, 
two  complementary  asymptotic  expansions  (useful  for  evaluating  the  Green 
function  sufficiently  far  away  from  the  singularity)  are  derived,  in  section 
6,  from  the  near-field  and  the  far-field  integral  representations  (5.11)  and 

(5.21) .  These  asymptotic  expansions  are  given  by  equations  (6.8)  and  (6.14). 
Comparison  of  these  complementary  expansions  then  yields  the  single  asymptotic 
expansion  (6.17).  In  section  7,  an  ascending  series  (useful  for  evaluating 
the  Green  function  in  the  vicinity  of  the  singularity)  is  derived  from  the 


near-field  integral  representation  (5.11).  This  ascending  series  is  given  by 
equations  (7.7),  (7.8),  and  (7.22).  In  section  8,  one-dimensional  Taylor 
series  expansions  are  obtained  from  the  near-field  and  the  Haskind  integral 
representations  (5.11)  and  (5.8c).  These  series  are  given  by  equations  (8.8) 
and  (8.13).  Finally,  expressions  for  the  gradient  of  the  Green  function  are 
given  in  section  9.  In  particular,  the  vertical  derivative,  Gz>  of  the  Green 
function  G  can  be  directly  expressed  in  terms  of  G,  as  shown  in  expression  (9.5), 
which  was  previously  given  by  Martin  [2],  and  is  obtained  here  in  a  different 
manner  (by  following  an  idea  used  by  Eggers  [4]  for  the  analogous  problem  of 
ship  wave  resistance). 
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The  basic  potential-flow  problem 

The  basic  potential-flow  problem  of  the  linearized  theory  of  flow 
about  a  body  in  regular  water  waves  investigated  in  the  present  study  is 
briefly  formulated  in  this  section.  A  sea  of  infinite  depth  and  lateral  ex¬ 
tent  is  assumed,  and  water  is  regarded  as  homogeneous  and  incompressible 
(with  density  p )  as  well  as  inviscid.  The  only  body  force  considered  is  that 
due  to  a  uniform  gravitational  field  (with  acceleration  g ).  Surface  ten¬ 
sion  and  free-surface  nonlinearities  are  neglected.  The  flow  is  irro- 
tational  and  thus  can  be  represented  by  a  velocity  potential  which  is 
a  function  of  the  Cartesian  coordinates  ^(X,Y,Z)  and  of  the  time  T,  i.e. 
<J>'(X,T).  The  mean  (undisturbed)  free  surface  of  the  sea  is  taken  as  the 
plane  Z  =  0,  with  the  Z  axis  positive  upwards. 

The  linearized  dynamic  free-surface  boundary  condition  takes  the  well- 
known  form 


gE'+*'+P'/p  =0  on  Z  =  0, 

where  E'(X,Y,T)  is  the  elevation  of  the  free  surface  above,  or  below,  its 
mean  level  Z  =  0,  4>'  =  84>'(X,Y,Z=0,T  )/3T,  and  P^(X,Y,T)  is  the  difference 

between  the  pressure  at  the  free  surface  and  the  atmospheric  pressure.  For 
most  problems  of  practical  interest,  the  pressure  at  the  free  surface  is 
constant  equal  to  the  atmospheric  pressure,  so  that  one  then  has  P'  =  0.  In 
the  presence  of  a  fluid  flux,  Q'(X,Y,T)  say,  across  the  free  surface,  the 
linearized  kinematic  tree-surface  boundary  condition  takes  the  form 

=  E'-Q'  on  Z  =  0, 

where  Q'<0  clearly  corresponds  to  fluid  being  sucked  away  across  the  free 
surface.  While  for  all  practical  problems  we  have  Q'  =  0,  it  will  be  impor¬ 
tant  to  allow  a  fluid  flux  across  the  free  surface  for  determining  the  free- 
surface  condition  verified  by  the  Green  function,  as  will  be  shown  in  the 
next  section.  Elimination  of  the  free-surface  elevation  E'  between  the  fore¬ 
going  dynamic  and  kinematic  free-surface  conditions  then  yields  the  free- 
surface  boundary  condition 


g$Z+$TT  =  ~P^/p-8Q'  on  Z  =  0, 


(1.1) 


which  thus  involves  the  velocity  potential  4>'  alone. 

In  the  present  study,  we  are  interested  in  flows  that  are  simple 
harmonic  in  time,  say  with  radiant  frequency  ^(period  2tt/u>).  As  it  is  well 
known,  and  is  discussed  for  instance  in  Stoker  [5],  such  free-surface 
gravity  flows  are  not  completely  (uniquely)  determined  unless  one  imposes  a 
"radiation  condition"  expressing  that  waves  at  a  sufficient  distance  away 
from  the  disturbance  (e.g.  a  body  in  the  present  problem)  which  created 
them  must  be  like  "outgoing"  progressing  waves,  i.e.  like  progressing  waves 
moving  away  from  the  wavt 3  source.  A  convenient  alternative  approach, 
employed  for  instance  in  Lighthill  [6] ,  to  the  use  of  such  a  "radiation  con¬ 
dition"  of  "outgoing  waves",  consists  in  defining  a  time-harmonic  flow  as 
the  limit-  as  the  small  positive  auxiliary  parameter  t  vanishes-  of  a  flow 
defined  by  a  velocity  potential  of  the  form 


<T(X,T)  =  Re$(X)  exp  [-ico(l+ie  )T] 


(1.2a) 


where  Re  represents  the  real  part.  The  free-surface  pressure  and  flux 
are  similarly  assumed  to  be  of  the  form 


P"(X,Y,T)  =  ReP(X,Y)exp[-iw(l+ie)T]  , 
Q'(X,Y,T)  =  ReQ(X,Y)exp[-iw(l+ie )T]  • 


(1.2b) 

(1.2c) 


In  this  alternative  approach,  one  then  is  faced  with  a  traditional  "initial- 
value  problem",  with  the  obvious  initial  conditions  $'  =  0  and  =  0  for 
T  =  -°°.  Use  of  expressions  (1.2a,b,c)  into  equation  (1.1)  then  yields  the 
following  free-surface  boundary  condition: 

g4-2-j2(l+ie)2$  =  iu)(l+ie  )P/p-gQ  on  Z  =  0,  (1.3) 

for  the  "spatial  component"  $(X)  of  the  actual  potential  $'($,T). 

It  will  be  convenient  to  define  adimensional  variables  in  terms  of 

1/co  as  reference  time  and  of  some  reference  length  L,  from  which  the  refer- 

2  2  2 

ence  velocity  <oL,  potential  <oL  ,  and  pressure  p<o  L  can  be  readily  formed. 
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We  thus  define  the  adimensional  variables 

t  =  U)T,  x  =  X/L,  <f>  =  <t>/u)L2,  p  =  P/pu)2L2,  q  =  Q/wL.  (1.4) 

In  terms  of  these  adimensional  variables,  the  free-surface  condition  (1.3) 
can  be  shown  to  become 

<J>  -f  (1+ie  )2<j>  =  if(l+ie)p-q  on  z  =  0,  (1.5) 

where  f  is  the  "frequency  parameter"  defined  as 

f  =  w2L/g.  (1.5a) 

The  "frequency  parameter"  f  can  obviously  be  made  equal  to  unity  by  select- 

2 

ing  the  reference  length  L  as  g/io  .  This  choice  of  reference  length 
essentially  corresponds  to  taking  the  length  of  the  water  waves  as  refer- 
ence  length  since  we  have  g/w  =  X/2tt  -  with  \  the  wave  length  of  plane  pro¬ 
gressive  waves  of  frequency  w  -  from  the  "dispersion  relation"  for  water  waves 
in  deep  water.  In  this  choice  of  reference  length,  the  size  of  the  body 
causing  the  waves  would  however  appear  to  vary  with  the  frequency  w  (the 
body  becoming  small  at  low  frequency  and  large  at  high  frequency  ).  An  al¬ 
ternative  (possibly  more  convenient  for  practical  purposes)  choice  is  to 
take  the  reference  length  L  as  a  length  characterizing  the  size  of  the  body, 
which  would  thus  remain  the  same  at  all  frequencies.  The  length  of  the 
waves,  however,  would  then  vary  with  the  frequency  (the  waves  being  long  at 
low  frequency  and  short  at  high  frequency). 

The  basic  potential-flow  problem  of  the  linearized  theory  of  flow 
about  a  body  in  regular  water  waves  can  now  be  briefly  stated.  As  it  is 
well  known,  this  problem  consists  of  solving  the  Laplace  equation 

V2<f>  =  0  in  (d),  (1.6a) 

subject  to  the  boundary  conditions  specified  below.  The  solution  domain 
(d)  in  equation  (1.6a)  is  the  domain  exterior  to  the  body  and  bounded  up¬ 
wards  by  the  mean  free  surface  (a)  say,  which  consists  in  the  whole  plane 


z  =  0  if  the  body  is  fully  submerged  or  in  the  portion  of  the  plane  z  =  0 
exterior  to  the  body  in  the  case  the  body  pierces  the  free  surface.  On  the 
mean  free  surface  (o),  the  free-surface  boundary  condition  (1.5)  must  be 
verified: 

<J>z~f  (1+ie)  2<j>  =  if(l+ie)p-q  on  (a),  (1.6b) 

where  in  fact  we  generally  have  p  =  0  =  q  for  the  problem  of  flow  about  a 
body.  The  potential  <J»  must  vanish  at  infinity;  specifically,  we  have  the 
condition 

p  =  0(l/|x|)  as  | x |  ->-oo  ,  (1.6c) 

expressing  that  p  vanishes  at  least  as  fast  as  1/ 1 x  |  as  |x|  -+  °°  .  Finally,  on 
the  body  surface,  (b)  say,  which  actually  consists  only  in  the  portion  of  the 
body  surface  located  below  the  plane  z  =  0  if  the  body  pierces  the  free  sur¬ 
face,  the  potential  <(>  must  verify  the  usual  "body  boundary  condition" 


<J>  given  on  (b)  , 


(1.6d) 


where  =  dp/dn  £  V<|>.n  is  the  derivative  of  p  in  the  direction  of  the  unit 
normal  n  to  (b) ,  which  is  taken  to  be  pointing  outwards  with  respect  to  the 
fluid,  i.e.  inside  the  body.  The  precise  form  taken  by  the  expression  for  <f>n 
on  (b)  in  particular  problems,  notably  in  the  usual  "radiation"  and  "diffrac¬ 
tion"  problems,  may  be  found  in  various  places  in  the  literature,  e.g.  in 
Newman  [7].  In  the  zero-frequency  limit  (f=0)  the  free-surface  boundary  con¬ 
dition  (1.6b)  becomes 


Pz  =  -q  on  (o)  ,  (1.7a) 

whereas  in  the  infinite-frequency  limit  (f=<*>)  it  becomes 


p  = -ip/ (1+ie)  on  (a). 


(1.7b) 
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2 .  The  Green  function 

A  well-known  technique  for  solving  a  potential-flow  problem  such  as 
the  one  defined  above  by  equations  (1.6  a,b,  c,  d)  ,  in  the  general  case  of  an 
arbitrarily-shaped  body,  consists  in  formulating  an  integral  equation  for  the 
potential  <£  based  on  the  use  of  a  Green  function  verifying  all  the  boundary 
conditions  of  the  problem  except  the  "body  condition",  which  is  to  be  verified 
by  means  of  the  integral  equation.  The  Green  function,  G(x;£,f,e)  say, 
appropriate  to  the  present  problem  then  is  the  solution  of  the  problem  de¬ 
fined  by  the  following  equations: 


7  G  =  5(x-£)  6(y-n)  6(z-0  in  z  <  0 


(2.1a) 


G  -f (1+ie)  G  =0  on  z  =  0, 
z 


G  =  0(l/r)  as  r  ->  °°  , 


(2.1b) 


(2.1c) 


where  <5(  )  is  the  usual  "Dirac  delta  function",  and  r  =  |x~C|  is  the  distance 
between  the  "field  point"  x  and  the  singularity"  £. 

A  particular  solution  of  the  Poisson  equation  (2.1a)  is  given  by 
4ttG  =  -1/r,  as  is  well  known  and  can  readily  be  verified.  The  general  solu¬ 
tion  of  equation  (2.1a)  can  thus  be  written  as 


4nG(x;£,f,e)  =  _l/r  +  H(x;£,f,e), 


(2.2) 


where  the  function  H  is  regular  harmonic  in  the  lower  half  space  z  <  0,  and 
evidently  is  to  be  determined  from  the  boundary  conditions  (2.1b ,c).  Indeed, 
use  of  expression  (2.2)  into  equations  (2.1a,b,c)  yields 


v  H  =  0  in  z  <  0, 


H  -f(l+ie)2H  =  [9  -f(l+ic)2] (1/r)  on  z  =  0, 
z  z 


H  =  0(1/ r)  as  r 


(2.3a) 

(2.3b) 

(2.3c) 


The  above  problem  can  be  solved  by  using  a  double  Fourier  transform 
with  respect  to  the  horizontal  coordinates  x  and  y.  The  double  Fourier  trans- 
form  of  the  function  H(x;£,f,e)  is  denoted  by  H  (a»B»z»£»f *e)  and  defined  as 
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It  may  be  seen  from  equation  (2.5)  that  the  first  term  on  the  right  side  of 

equation  (2.7b)  is  equal  to  the  double  Fourier  transform  (1/rO  of  1/r*, 

2  2  21/2 

where  r '  is  defined  as  r'  =  (x'  +y'  +z'  )  with  x'  =  x-£,  y'  =  y-n ,  and 
z ’  =  z+c ,  so  that  x'(x',y',z')  is  the  vector  joining  the  mirror  image  of  the 
"singularity"  £  with  respect  to  the  free  surface  z  =  0  to  the  "field  point" 
x,  and  r'  is  the  distance  between  these  two  points. 

The  function  H(x;£,f,e)  may  now  be  obtained  by  taking  the  inverse  double 

Jcic  -v 

Fourier  transform  of  the  function  H  (a ,8 ,z ;£ , f ,e) ,  namely 

.  i  f  r  -i(ax+By)  **  ^ 

H(x,y  ,z ;  f  ,e)  =  dB  I  da  e  H  (a,  8  ,z  ;£  ,f ,  e)  .  (2.8) 

'—CO  /— OO 


By  using  expressions  (2.7a,b)  for  H  into  equation  (2.8),  and  by  using  the 
resulting  expression  for  H  into  equation  (2.2),  we  can  then  obtain  the  follow¬ 
ing  alternative  expressions  for  the  Green  function  G(x;f,f,e): 


4nG  =  -  — 
r 


1  _  f (1+ie) 2 
r '  7r 


r 


L 


d8 


00  '—Oo 


f°°  kz  ^-i(ax'+8y  ') 

da  e 


k[k-f(l+ie)  ] 


(2.9a) 


4ttG  = 


i  +  f 

r  r  it  J 


d6 


kz^-i(ax^+8y') 

da  e  _ 

/_oo 


k-f (1+ie) 


(2.9b) 


The  Green  function  G(x;£,f,e)  obviously  is  axisymnetric  about  the 

vertical  axis  x=£,  y=n ,  so  that  we  may  take  y'  =  y-n  as  zero  and  replace 
2  2  1/2 

x'  =  x-£  by  p  =  (x'  +y'  )  in  expressions  (2.9a,b).  Expression  (2.9b)  then 
becomes 


4ttG 


..  22. 1/2  . 
z  (a  +8  )  -ipa 

da  e _ 

(a2+B2)1/2-f(l+ie) 


2 


(2.10a) 


A  more  usual  alternative  expression  Is  that  which  can  be  obtained  by  performing 
"M‘  ll,,t  ,on  y  ’  -  0  and  x  ‘  =  ()  in  expression  (2.9n),  followed  by  a 


transformation  from  the  Cartesian  Fourier  variables  a  and  £  to  the  polar 
2  2  1/2 

variables  k  =  (a  +3  )  and  0,  specifically  by  performing  the  change  of 
variable  a  =  kcosB  and  0  =  ksin0.  The  resulting  classical  expression  is 

,  ,  oof"  (z'-ipcosO)k 

4uG  =  -  —  -  f(l+ie)z  -  f  do  |  dk  e  (2.10b) 

r  r  11  J  J  5  • 

0  u  k-f (1+ie)2 

Expressions  (2.10a,b)  show  that  G(x;t»f.e)  actually  is  a  function  of  only 

2  2  1/2  2 

three  space  variables,  namely  p  =  (x"  +y'  )  ,  z'  =  z+4,  and  (z-4)  which 

2  2  1/2 

occurs  in  r  =  [p  +(z-4)  J 

As  it  is  self-evident  from  equations  (2.1a,b,c),  the  physical  signifi- 
cance  of  the  Green  function  G(x;£,f,e)  is  that  ReG(x;£ ,f ,e)exp(e-i) t  is  the 
linearized  velocity  potential,  at  the  "field  point"  x(x,y,z<0)  and  at  the 
time  t,  of  the  flow  caused  by  a  submerged  pulsating  source  of  strength 
Reexp(e-i) t ,  -»  _<  t  <_  t,  located  at  point  t(S»n»C<0).  This  well-known 
physical  interpretation  becomes  ambiguous  however  in  the  limiting  case  4  =  0, 
since  the  source  then  is  obviously  no  longer  fully  submerged.  A  natural  com¬ 
plementary  interpretation  for  this  limiting  case  is  to  assume  that  the  out¬ 
flow  produced  at  point  (£,n,4=0)  now  stems  from  a  flux  Req(x,y)exp(e-i)t , 
with  q(x,y)  =  6(x-£) 6 (y— n ) ,  across  the  mean  free  surface  z  =  0.  Equations 
(1.6a,b,c)  then  suggest  that  the  "limit  Green  function"  G£(x;£,n,f,e)  = 

G(x; £,n ,4=0 , f  ,e)  must  verify  the  following  equations 

V2G^  =  0  in  z  <  0,  (2.11a) 

G^z~f (1+if ) 2G^  =  -6(x-£)6(y-n)  on  z  =  0 ,  (2.11b) 

G^  =  0(l/r)  as  r  -*  °°  .  (2.11c) 

A  mathematical  demonstration  of  the  above  physically-motivated  equa¬ 
tions  can  readily  be  provided  by  verifying  that  the  solution  G^(x;£,n,f ,e)  of 
the  problem  defined  by  equations  (2.11a,b,c)  actually  is  identical  to  the  "limit 
Green  function"  obtained  by  replacing  4  by  zero  in  the  previously-derived 
solution  G(x;ff,f,e)  of  the  problem  defined  by  equations  (2.1a,b,c).  Problem 
(2.11)  may  be  solved  in  the  same  manner  as  was  used  previously  for  solving 
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problem  (2.3),  namely  by  using  a  double  Fourier  transform  with  respect  to 
the  horizontal  coordinates  x  and  y.  We  may  then  obtain  the  "Fourier-trans- 
formed  problem: 

d2G^*/dz2-k2G^*  =  0  in  z  <  0, 

dG^  /dz-f  (1+ie)  G^  =  -exp  [ i )  ]/2tt  on  z  =  0, 

G^  -*•  0  as  z  -y  -  , 

where  G^  is  the  double-Fourier  transform  of  G^ ,  as  is  defined  by  formula 
(2.4).  The  solution  of  the  above  problem  is  given  by 

**  2 
G^  =  -  exp  [kz+i(a?+6n)  ] / 2tt [k— f  (1+ie)  ] 

By  taking  the  inverse  double  Fourier  transform,  as  is  given  by  formula  (2.8), 
we  can  finally  obtain 


G^(x;C,n,f,e) 


•o° 

dS 

—00 


kz-i(ax'+By ') 
da  e _ 

k-f (1+ie)2 


(2.12) 


which  can  readily  be  verified  to  be  identical  to  the  expression  obtained  by 
replacing  4  by  zero  in  formula  (2.9b). 

Conversely,  it  may  be  shown  that  the  "limit  Green  function"  G^  given 
by  expression  (2.12)  does  in  fact  verify  equations  (2.11a,b,c).  Verification 
of  equations  (2.11a)  and  2.11c)  can  easily  be  checked.  As  for  the  free- 
surface  condition  (2.11b),  we  have 


G£z-f(l+ic)2G£ 


—00 


-iax'  1 
e  da)  (^ 


r°°  -iBy' 

e  dB) 


on  z  =  0, 


from  which  we  may  obtain 

G^-f  (1+i'.)  2G^  =  -6(x')6(y')  =  — cS  ( x— f;)  <S(y-n)  on  z  =  0 


by  virtue  of  the  relations 


1  = 


6(x)dx,  6 (x)  = 


-iax 
e  da, 


expressing  the  (well-known)  fact  that  6(x)  and  1  are  Fourier  transforms. 

It  may  thus  be  seen,  in  summary,  that  the  Green  function  G(x;£,f,c) 
of  the  theory  of  flow  about  a  body  in  regular  waves  (where  the  limit  e  =  +0 
is  ultimately  implied)  verifies  the  following  equations 

(2.13a) 
(2.13b) 
(2.13c) 

(2.14a) 
(2.14b) 
(2.14c) 

As  it  was  noted  previously,  the  Green  function  only  depends  on  the 
2  2  2 

three  space  variables  (x-£)  +(y-n)  ,  (z~0  ,  and  (z+c),  so  that  this  function 
is  invariant  under  the  substitution  x  it.  Physically,  the  velocity  poten- 
tial  ReG(x;C,f ,e)exp(E-i)t  of  the  flow  created  at  point  x(x,y,z<0)  by  an  out¬ 
flow  of  strength  Reexp(e-i)t  at  point  £(£,n,c< 0) ,  stemming  from  a  submerged 
source  if  ?<0  or  a  free-surface  flux  if  ?  =  0,  is  identical  to  the  potential 
ReG(£ ;x , f , e) exp(e-i) t  of  the  flow  created  at  point  £  by  an  outflow  Reexp(e-i)t 
at  point  x,  s  emming  from  a  source  if  z<0  or  a  free-surface  flux  if  z  =  0. 

It  then  follows  that  equations  (2.13)  and  (2.14)  are  also  verified  by  the 
function  G(£;x, f,e). 

In  the  zero-  and  infinite-  frequency  limits  (f=0,°°),  the  Green  func¬ 
tion  G  becomes 


v  G  =  6(x-C)6(y-n)S(z-c)  in  z  <  0, 
Gz-f(l+ie)2G  =0  on  z  =  0, 

G  =  0(1/ r)  as  r  +  “  , 


if  C  <  0 


v  G  =  0  in  z  <  0, 

Gz-f(l+i£)2G  =  -6 (x-O 6 (y-n)  on  z  =  0, 
G  =  0(l/r)  as  r  -*  «  . 


if  4  =  0 


4-G°  =  -1/r-l/r' 


4HG"  =  -1/r+l/r ' 


(2. 15a,b) 


respectively,  as  may  be  obtained  from  formulas  (2.10a)  and  (2.10b),  respectively. 
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The  zero-frequency  Green  function  G^(£;x)  verifies  the  equations 


V2G°  =  <5(x-£) 6(y-n) 6(z-t)  in  z  <  0, 


G  =0  on  z  =  0, 
z 


2  0 

V  G  =0  in  z  <  0, 

„0 


if  t  <  0, 


>  if  5  =  0, 


(2.16) 


-6(x-C)  6(y-ri)  on  z  =  0  , 

€t 

whereas  the  infinite-frequency  Green  function  G  (£;x)  verifies  the  equations 


V  G  =  6(x-C)  <5 (y— ri)  S(z-t)  in  z  <  0, 
G  =0  on  z  =  0, 

G"  =  0  if  c  =  0. 


V  if  C  <  0, 


(2.17) 
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3.  Basic  Integral  identities 

In  this  section,  basic  integral  identities  for  the  velocity  potential 
4(x)  are  obtained  by  applying  a  classical  Green  identity  to  the  function  4(x) 
and  the  Green  function  G(C;x,f,e)  defined  in  the  previous  section.  The  above- 
mentioned  classical  Green  Identity  is 


(<JiV2G-GV2({>)dv 

d 


(<(>G  -G<t>  Jdxdy  +  (4G  -G4„)da, 

>o  Z  >b  n  n 


(3.1) 


where  the  fact  that  <}>  =  0(l/r)  and  G  »  0(l/r)  as  r  -*■  »  was  used  to  discard  the 

3 

integral  of  4G  — G<}>  =  0(l/r  )  over  a  large  surrounding  half  sphere  of  radius  r 

n  2n 

(surface  area  ^  r  )  .  In  formula  (3.1),  as  indeed  hereafter  in  this  study, 

the  following  notation  is  used:  4>  =  4(x)  ,  G  =  G(it;x,f  ,e) ,  V  =  Ox»3y»^z^* 

G  =  9G/9z,  G  =  9G/9n  =  VG-n,  and  likewise  41  =  94/9n  =  74 -n,  where  ft  =  n(x) 
z  n  n 

is  the  unit  normal  vector  at  point  x  of  the  body  surface  (b)  pointing  outside 
the  solution  domain  (d)  as  was  defined  in  section  1;  furthermore,  dv  =  dv(x) 
and  da  =  da(x)  represent  the  differential  elements  of  volume  and  area  at 
point  x  of  (d)  and  (b)  respectively,  while  dxdy  evidently  is  the  differential 
element  of  area  of  the  mean  free  surface  (a) . 

2  2 

By  expressing  the  term  <j>G  — G<J>  in  the  form  4[G  -f(l+ie)  G ]  — G [ <+>  -f(l+ie)  4], 

z  z  z  z 

equation  (3.1)  can  be  written  as 


(  2  I  2  (  2  (  2  f 

4>V  Gdv  -  4>[Gz-f(l+ie)  G]dxdy  =  GV  <j>dv-  G[4z-f(l+ie)  4]dxdy-l  (G4  -<t>G  )da.  (3.2) 

■'d  ^  c  J  d  '  o  •'b 


Let  us  now  express  4  in  the  two  integrals  on  the  left  side  of  equation  (3.2) 
in  the  form  4*+<t>~4*  ,  where  $  =  4(x)  as  was  defined  previously,  and  4.  is 

-V  ->-  * 

defined  as  4(C).  i-e.  4*  -  4(C)-  We  may  then  obtain 
(of  0 

4V^Gdv-|  4 [ G  -f (1+ie)  G]dxdy  =  C4*+<T  ,  (3.3) 

Jd  '  cr  2 

where  C  and  C'  are  defined  as 


C  = 


2 

V  Gdv- 
d 


[Gz~f (1+ie) 2G] dxdy , 


(3.4) 


0 
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>d 


(<Hf>*)7' Gdv  - 


f(l+ie)  G]dxdy. 


It  follows  from  equations  (2.13a,b)  and  (2.14a,b),  and  the  fact  that  =  0 

when  the  points  x  and  £  coincide,  that  we  have  C'  r  0  for  any  point  £  in  the 
lower  half  space  T,<_ 0.  By  using  equation  (3.3),  with  C  '  5  0,  into  equation 
(3.2),  we  can  then  obtain 

=  |  GV2 $dv  -  G [ -f  (1+ie)  2<j>]dxdy  -  I  (Gd>  -4>G  )da,  (3.5) 

1  d  >o  Z  •'b  n  n 


where  C  is  given  by  formula  (3.4). 

By  using  equations  (2.13a,b)  and  (2.14a,b)  in  expression  (3.4)  for 
C,  it  can  be  seen  that  we  have  C  =  1  if  the  point  £  is  strictly  outside  the 
body  surface  (b) ,  either  in  (d)  or  on  (o) ,  while  we  have  C  =  0  if  £  is  strictly 
inside  (b) .  We  thus  have 


<t>*  =  |  GV2<fidv-j  G[<f>z- f(l+ie)2<p]dxdy  ~ 

|b(GV-,Gn>da 

(3.6a) 

for  £  in  (d)+(o)-(b)-(c) ,  where  (c)  is 

the  intersection  curve  of  (b)  with  the 

plane  z=0  [in  the  case  (b)  does  in  fact  intersect  the  plane  z=0],  and 

f  o  (  2  ( 

0  «  !  G5/  4>dv  -  G [ <J>  -f (1+ie)  <f>]dxdy  -  1 

^d  >c  2  j 

(G<f>  -4>G  )  da 
,  n  n 

b 

(3.6b) 

for  £  in  (di)+(ai)-(b)-(c) ,  where(di)  and  (cO  represent  the  domain  and  the 
portion  of  the  plane  z=0,  respectively,  which  are  strictly  inside  the  body 
surface  (b)+(c) ,  so  that  (d)+(dj)+(b)  and  (a)+(ai)+(c)  are  the  whole  lower 
half  space  z< 0  and  the  entire  plane  z=0,  respectively.  It  can  also  easily 
be  seen  from  equations  (2 .13a,b) , (2. 14a, b) ,  and  (3.4)  that  we  have  C=l/2  if 
the  point  f  is  right  on  the  body  surface  (b)  or  its  intersection  (c)  with 
the  plane  z=0,  at  least  for  points  %  where  (b)+(c)  is  smooth;  more  generally 
the  value  of  4irC{or  2ttC >  at  a  point  £  of  (b)  (or  (c)>  is  equal  to  the  angle 
under  which  (d)  (or  ( o) }  is  viewed  from  the  point  £•  We  then  have 


« 
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I - 

i  ■■!** 

i _ 


GV2<j>dv  - 

'd 


f  o 

I  G [ 4>z — f ( 1-Hie )  <f>]dxdy  - 

'  a 


(G<p  — 4>G  ) da 
n  n 


(3.6c) 


for  t  exactly  on  (b)+(c).  Equations  (3.6a,b,c)  are  quite  well  known  of 

course,  although  the  particular  case  of  these  equations  corresponding  to 
2  2 

V  if,  =  0  =  <J>z-f(l+ie)  <(>  is  usually  given  in  the  literature,  and  the  tradi¬ 
tional  derivation  of  these  equations  from  the  Green  identity  (3.1)  is  some¬ 
what  different  from  that  shown  above. 

As  it  is  explicitly  indicated  in  equation  (3.6a,b,c),  the  value  of 
the  constant  C  on  the  left  side  of  equation  (3.5)  is  discontinuous  across 
the  body  surface,  C  being  equal  to  1  outside  the  body  and  to  0  inside.  This 
discontinuity  in  the  value  of  C  on  the  left  side  of  equation  (3.5)  evidently 
is  accompanied  by  a  corresponding  discontinuity  on  the  right  side  of  the 
equation.  Specifically,  the  latter  discontinuity  stems  from  the  integral 
da  representing  the  potential  induced  by  a  surface  distribution  of 
normal  doublets  of  strength  <)>  over  (b) ,  as  is  well  known.  An  integral  iden¬ 
tity  valid  for  any  point  t — either  outside,  inside,  or  right  on  the  body 
surface — can  be  obtained  by  eliminating  the  discontinuity  in  the  value  of 
C  in  equation  (3.5).  This  can  readily  be  achieved  by  adding  the  term 
on  both  the  left  and  right  sides  of  equation  (3.5),  with  given  by 

[Gz-f(l+ie)2G]dxdy. 


Equation  (3.5)  then  becomes 


1 

*■  ✓  j 


GV2<J>dv  -  Gt^-f (l+ie)2<t>]dxdy  -  j  (G<t>n-<t>Gn)da  + 

d  o  -"b 


+<(.*{ 


V2Gdv 


[Gz-f(l+ie)2G]dxdy} 


(3.7) 


where  I  is  defined  as 


I  = 


2  f  2 

V  Gdv  -  I  [G  -ffl+ie)  Gjdxdy  . 

^d+d^  -*0+0^ 
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It  can  be  seen  from  equations  (2.13a,b)  and  (2.14a,b)  that  we  have  1=1  for 
any  point  £  in  the  lower  half  space  £<0 .  By  using  the  divergence  theorem 


(3.8) 


in  equation  (3.7),  where  1  is  taken  as  1,  we  may  finally  obtain 


f  2  f 

(l-w^H*  =  GV  <j>dv  - 

Jd 

G [<J>  -f  (1+it)  2<J>]dxdy  -  [G<f>  — Cq>— . ) G  ]da 

_  z  ,  n  *  n 

a  'b 

where  the  function  =  w(£;f,E,cK)  is  the  "waterplane  integral" 

w(t;f,e,oi)  =  f(l+ie)2 

' 

G(£;x,f ,e)dxdy 

'  0  . 

1 

• 

,  (3.9) 


(3.9a) 


An  interesting  alternative  form  of  the  integral  identity  (3.9)  involving  a 
"modified  Green  function  "  G  can  be  obtained  by  dividing  both  sides  of  this 
identity  by  (1-w^) .  This  yields 


=  (  &v24>dv  -  (  &[<t>z-f(l+i£)2()>]dxdy  - 
'  d  '  a 

where  the  "modified  Green  function"  C!  is 

given  by  £  =  G/(l- 

d(l;^,f,e,a.)  =  G(£;J,f ,e)/[  1-f  (1+ie)  2  j 

G(£;x,f ,e) dxdy ] 

°i 

(3.10) 


that  is 


(3.10a) 


In  the  case  of  a  fully-submerged  body,  we  evidently  have  (o^)  =0,  so  that 
the  "waterplane  integral"  w  then  vanishes  and  the  "modified  Green  function" 

G  becomes  identical  to  the  usual  Green  function  G,  that  is  we  have  w  =  0  and 

%  =  G. 

The  integral  identity  (3.9)  is  valid  for  any  point  £  either  outside, 
inside,  or  exactly  on  the  body  surface.  This  new  identity  thus  is  essenti¬ 
ally  equivalent  to  the  set  of  the  three  classical  identities  (3.6a,b,c), 


which  are  valid  exclusively  for  £  outside,  inside,  or  on  the  body  surface, 
respectively.  As  a  matter  of  fact,  these  three  identities  can  readily  be  ob¬ 
tained  from  the  identity  (3.9)  by  verifying  that  we  have  w  -  1  G  da  =  0,1,  or 
1/2  for  £  outside,  inside,  or  on  the  body  surface,  respectively,  as  is  shown 
below.  From  equation  (3.9a) ,  we  have 


f  2 

w .  -  G  da  =  f(l+ie)  Gdxdy  -  G  da  = 
*  /,  n  ,  n 

•  h  J  n  h 


G  dxdy  -  G  da 

o .  '  b 

i 


(Gz-f(l+ie)  Gldxdy 


(  j  f  2 

V  Gdv  -  [Gz-f(l+iE)^G]dxdy, 

d.  '  a . 


where  the  divergence  theorem  (3.8)  was  used;  the  above-stated  result  then 
follows  from  equations  (2.13a,b)  and  (2.14a,b). 

Although  the  present  study  is  primarily  concerned  with  the  "exterior 
potential-flow  problem",  i.e.  the  problem  of  potential  flow  about  a  body,  it 
may  be  of  some  interest  to  list  here  the  integral  identities  corresponding 
to  equations  (3.6a,b,c)  and  equations  (3.9)  and  (3.10)  for  the  "interior 
potential",  <)>*  say,  defined  in  the  interior  domain  (d^)  bounded  by 
(b)+(a^)+(c) .  The  integral  identities  corresponding  to  equations  (3.6a,b,c) 
can  be  shown  to  take  the  forms 


<j>*  ■  GV^ifi^dv  -  Gf^*  -  f (1+ie) ^^Idxdy  + 
>A  Jr,  z 


"7;  1 

(G**-*  G)da  | 
n  n 


(3.11a) 


for  t  in  (d^)+(o^)-(b)-(c) ,  that  is  inside  the  body. 


f  2  i 

0  *  GV  $  dv 

J  d. 


-f  gi*7 

*  rr 


f  (1+ie)  dxdy  + 


Ti  1 

(G*  -$xG)da  I 
n  n 


(3.11b) 


for  £  in  (d)+(o)-(b)-(c) ,  i.e.  outside  the  body,  and 


GV  d>  dv  - 


G[(|>  ~f(l+ie)  <j>  ]dxdy  +  (G<J)  -<|>  G  )da 

>n  Z  Jh  n  n 


(3.11c) 


for  £  exactly  on  (b)+(c) . 

The  integral  identity  corresponding  to  equation  (3.9)  takes  the  form 

=  (  G72<(.1dv  -  (  G[<(>^-f (l+ic)2(j)1]dxdy  +  [ G^1  - ( ch 1 1) G  ]da  I  ,  (3.12) 

Jh  Jn  Z  Jh  n  n  I 


where  the  function  w.  =  w(£;f,e,o.)  is  the  "waterplane  integral"  defined  by 

X  1 

formula  (3.9a).  The  integral  relation  (3.12),  like  equation  (3.9),  is  valid 
for  any  point  t,  either  inside,  outside  or  exactly  on  the  body  surface,  and 
thus  is  equivalent  to  the  set  of  the  three  classical  identities  (3.11a,b,c). 
An  alternative  form  of  the  integral  identity  (3.12)  involving  a  "modified 
Green  function"  G*  is 


G^V^^dv  -  |  Gi  [ <f> f  (1+ie)  2<f>'i']dxdy 

a. 

l 


i 

+  [Gi^^-(d>±-<(>i)Gf  ]da  j  ,  (3.13 


where  the  function  G  is  given  by  G  =  G/w.  ,  that  is 


G1(I;x,f ,e ,oi)  =  G(|;x,f ,e)/f (1+ie)  G(|;x,f ,e)dxdy 


(3.13a) 


The  case  of  a  fully-submerged  body  evidently  is  a  degenerate  case  for  which 
the  boundary  condition  at  the  free  surface  becomes  meaningless. 

It  is  interesting  to  note  that  whereas  the  classical  identities 
(3.6a,b,c)  and  (3.11a,b,c)  for  the  exterior  and  interior  potential-flow 
problems,  respectively,  are  identical  (provided  only  one  properly  interchanges 
the  terms  "exterior"  and  "interior"),  the  relations  (3.9)  and  (3.12),  or  the 
equivalent  alternative  relations  (3.10)  and  (3.13),  clearly  are  not  identi¬ 
cal.  Specifically,  the  factor  (1-w^)  in  relation  (3.9)  for  the  exterior 
problem  becomes  w^  in  relation  (3.12)  for  the  interior  problem,  and  the 
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equivalent  relations  (3.10)  and  (3.13),  although  identical  in  form,  involve 
different  "modified  Green  functions",  namely  the  "modified  exterior  Green 
function"  G  H  G/(l-w^)  and  the  "modified  interior  Green  function"  G*  =  G/w^. 

If  we  add  the  integral  relations  (3.6a)  and  (3.11b),  we  may  obtain 
the  relation 


GV^<j)dv  + 
>  d 


2  i 

GV  <(>  dv  - 

d. 

i 


G[4>7~f  (1+ie)  %]dxdy 

o 


G  [  f  (1+if)  ]dxdy 
i 


+ 


[G(^-4>„)-(4>1-4>)G.  ]da. 
i  n  n  n 

b 


Addition  of  equations  (3.6b)  and  (3.11a)  yields  the  same  relation,  with  <t>*  on 
the  left  side  merely  replaced  by  4>^.  We  therefore  have  the  relation 


GV  <t>dv  - 


L  _  ^.<0_ 


z=0 


G [ 4>z— f (1+ie)  <(>]dxdy  + 


[G(*n"')>n)"('!'1"t)Gnlda  I  »  (3,14) 


where  41  on  the  left  side  and  in  the  first  two  integrals  on  the  right  side 
clearly  corresponds  to  if1  or  for  points  outside  or  inside,  respectively, 
the  body  surface  (b)+(c) .  The  integral  relation  (3.14)  expresses  the  poten¬ 
tial  <f>(£)  in  the  lower-half  space  C.<0  in  terms  of  a  volume  distribution  of 

2 

sources  (with  density  V  4>)  ,  a  free-surface  flux  (or  pressure)  distribution 

2 

[with  flux  q=f(l+ic)  4>-<j>  ],  and  surface  distributions  of  sources  (with  den- 

sity  4>n— d>n)  and  of  normal  dipoles  (with  density  4>-<f>  )over  the  surface  (b). 

Two  classical  results  in  potential  theory  immediately  follow  from  equation 

(3.14):  namely,  (i)  a  distribution  of  normal  dipoles  (say  of  strength  6)  on 

a  surface  (say  S)  generates  a  potential  4>  =  f  <5G  da  whose  value  is  discon- 

JS  n  e  i 

tinuous  across  the  surface  S  (specifically,  we  have  41  -<f>  =6,  where  the 
superscripts  e  and  i  refer  to  the  "exterior"  and  "interior"  sides  of  S,  re¬ 
spectively;  the  "interior"  side  being  that  into  which  the  unit  normal  vector 
n  is  drawn) ,  and  (ii)  a  distribution  of  sources  (strength  a)  over  a  surface 

S  generates  a  potential  <f>  =  J  Goda  whose  normal  derivative  4  is  discontin- 

6 

uous  across  S  (specifically)  we  have  4>n-4|n  =  o)  .  It  may  finally  be  worth¬ 
while  to  verify  that  equation  (3.14)  can  (naturally)  also  be  obtained  by 
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adding  equations  (3.9)  and  (3.12).  Indeed,  this  yields 


Gv  <j>dv  - 


z<0 


G [ cf>  -f(l+ie)  <()]dxdv  + 
i  z 


[G(<f,n  ' <’n)"(*1~<t’)Gnlda’ 


where  C  is  defined  as 


C. 

i 


f  f  1  \  _  f(d.)+(o.)-(b)-(c) 

G  da  =  |  )  for  (  in  i 

bn  L  0  J  [(d)+(a)-(b)-(c) 


as  was  shown  previously.  It  then  follows  that  the  expression  (1-C  ) <(>  +C  <t>* 

,  i  ■+  11 

is  identical  to  <f>A  or  <|>^  for  £  outside  or  inside,  respectively,  the  body 

surface  (b)+(c) ,  and  the  above  relation  thus  is  identical  with  relation  (3.14). 

In  the  zero-frequency  limit  (f=0),  the  integral  identity  correspond¬ 
ing  to  identity  (3.9)  takes  the  form 


G°V2<f>°dv  - 


GG<(>Gdxdy 


[G°d»°-(<J>0-4.°)G°]da 
n  *  n 


(3.15) 


where  G^  is  the  zero-frequency  Green  function  given  by  formula  (2.15a),  and 
is  the  zero-frequency  potential,  which  verifies  a  free-surface  condition  like 
equation  (1.7a).  The  zero-frequency  integral  identity  (3.15)  may  readily  be 
obtained  as  the  zero-frequency  limit  of  the  integral  identity  (3.9),  or  it 
may  be  derived  directly  by  using  equations  (2.16)  in  the  Green  identity  (3.1). 
In  the  infinite-frequency  limit  (f=°°),  on  the  other  hand,  the  integral  identity 
corresponding  to  identities  (3.9)  and  (3.15)  takes  the  form 


(3.16) 


OO  OO 
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oo  2  00 

G  V  dv  + 
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<f>  G  dxdy  - 
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r 

00  00  m  00.  00 

[G  «f>  -(♦  -dOGlda 
,  n  *  n 

b 

as  may  be  obtained  by  using  equations  (2.17)  in  the  Green  identity  (3.1). 

oo 

In  equation  (3.16),  G  is  the  infinite-frequency  Green  function  given  by 

OO 

formula  (2.15b),  <f>  is  the  infinite-frequency  potential,  which  verifies  a 

00  oc  H*- 

free-surface  condition  like  equation  (1.7b),  and  the  function  =  w  (£;a^) 
is  defined  by 


It  may  easily  be  seen  that  all  the  terms  in  equation  (3.16)  are  continuous 
at  the  body  surface  (b)  if  c<0.  It  can  also  be  seen  that,  in  the  limit  £=-0, 
the  first  and  last  integrals  on  the  right  side  of  equation  (3.16)  vanish, 

oo 

while  the  function  w*  and  the  second  integral  on  the  right  side  (i.e.  the 

two  free-surface  integrals  involving  G  )  are  discontinuous  at  the  intersec- 

z 

tion  curve  (c)  of  the  body  surface  (b)  and  the  mean  free-surface  plane  £=0; 
specifically,  the  terra  (1-w^  )<P *  on  the  left  side  and  the  integral  f  ^^G^dxdy 
may  be  shown  to  be  equal  to  0,<|>”/2,  and  for  f,(£,n,C=-0)  on  (oi),(c) 
[wherever  (c)  is  smooth],  and  (a),  respectively. 
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In  the  particular  case  when  we  do  not  have  a  body  (b) ,  equations  (3.6b) 
and  (3.6c)  evidently  become  meaningless,  while  both  equations  (3.6a)  and 
(3.9)  can  readily  be  seen  to  become 

<j>.  =  GV^<j>dv  -  G[<(>  -f(l+ie)2<t>]dxdy 

h<0  J2=o  z 


This  formula  provides  an  explicit  solution  to  the  problem  of  determining  the 

2  2 

velocity  potential  4>  when  V  <j>  and  4>  -f(l+ie)  (j>  are  given  in  the  lower-half 

z 

space  z<0  and  on  the  mean  free  surface  z=0,  respectively.  A  classical  prob¬ 
lem  in  this  category  is  that  of  determining  the  potential  induced  by  a  given 
distribution  of  pressure  p  (and  possibly  also  of  flux  q)  at  the  free  surface 
The  solution  of  this  problem,  which  is  stated  in  differential  form  by  equa¬ 
tions  (1.6a,b,c),  is  then  given  by 


G(€;x,y,z=0,f ,e) [if (l+ie)p(x,y) -q(x,y) ]dxdy 


This  solution  is  well  known  of  course,  and  can  indeed  be  obtained  by  using 
other  approaches  than  that  adopted  here.  For  instance,  equations  (1.6a,b,c) 
can  be  solved  directly  by  means  of  a  double  Fourier  transformation  with  re¬ 
spect  to  the  horizontal  variables  x  and  y.  Another  often-used  approach  for 
determining  the  potential  induced  by  a  free-surface  pressure  distribution  is 
based  on  the  use  of  a  special  Green  function,  G  say,  corresponding  to  the 
potential  induced  by  a  pressure  impulse  p(x,y)  =  6(x-£;)  S(y-n)  at  the  free 
surface.  This  Green  function  G  in  fact  is  essentially  identical  to  the  "limit 
Green  function"  G^fmore  precisely,  we  have  G  =  -if(l+ie)G£]  defined  by  equa¬ 
tions  (2.11a,b,c)  and  given  by  formula  (2.12).  The  "free-surface  pressure 
impulse"  Green  function  G  thus  is  merely  a  particular  case  (corresponding  to 
z=0)  of  the  Green  function  G  defined  in  this  study,  so  that  introduction  of 
this  special  Green  function  G  is  obviously  not  necessary. 

A  particular  problem  that  is  more  difficult,  and  also  more  important 
from  the  point  of  view  of  practical  applications,  than  the  above-discussed 
problem  of  flow  induced  by  a  free-surface  pressure  dltribution  is  the  problem 
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2 

of  flow  about  an  arbitrary  body  in  regular  waves,  for  which  we  have  V  4>  =  0 

2 

in  the  mean  flow  domain  (d)  and  p  =  0  =  q,  and  hence  <t>z~f(l+ie)  4>  =  0,  on  the 
mean  free  surface  (o)  in  the  usual  linearized  approximation.  The  integral 
identities  given  in  the  previous  section  evidently  do  not  provide  an  explicit 
solution  for  this  problem  as  they  did  in  the  previous  case  of  flow  induced 
by  a  free-surface  pressure  distribution  in  the  absence  of  a  body.  However, 
identities  (3.6c)  and  (3.9)  provide  integral  equations  for  determining  the 
potential  ^  =  $(£)  on  the  body  surface  (b)  +  (c) .  These  integral  equations 
take  the  form 


(4.1) 


(4.2) 


corresponding  to  equations  (3.6c)  and  (3.9),  respectively .  An  alternative 
form,  corresponding  to  equation  (3.10),  of  the  integral  equation  (4.2)  is 


(4.2a) 


A  notable  feature  of  the  new  integral  equation  (4.2),  or  of  the 
equivalent  equation  (4.2a),  is  that  the  last  integrals  on  the  right  sides  of 
these  equations  are  continuous  functions  of  X  (since  the  dipole  strength  Q-Q* 
vanishes  as  the  "integration  point"  x  and  the  "field  point"  X  coincide) , 
whereas  the  corresponding  integral  in  the  classical  integral  equation  (4.1) 
is  discontinuous  at  the  body  surface  (b)  +  (c) ,  as  was  discussed  in  the  pre¬ 
vious  section  and  is  indeed  well  known.  As  a  matter  of  fact,  the  factor  1/2 
on  the  left  side  of  equation  (4.1)  is  correct  only  at  points  X  where  (b)  +  (c) 
is  smooth,  that  is,  has  a  tangent  plane,  as  was  also  discussed  in  the  previous 
section.  The  classical  integral  equation  (4.1)  requires  evaluation  of  a  dis¬ 
continuous  function  (the  function  defined  by  the  last  integral  on  the  right 
side)  precisely  on  the  surface  of  discontinuity,  namely  (b)  +  (c) .  This 
awkward  and  difficult  problem,  specially  from  the  viewpoint  of  numerical 
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calculations,  is  avoided  in  the  integral  equations  (4.2)  or  (4.2a).  In  fact, 

*\/ 

the  integrands  ( <j>— <J>  )G  or  (ij>-<j>.)G  can  be  shown  to  be  non-singular  (that  is, 
to  remain  finite)  for  x  =  £,  which  evidently  is  advantageous  for  purposes  of 
numerical  calculations. 

For  simplicity,  let  the  first  integral  on  the  right  side  of  the 
integral  equation  (4.2)  be  denoted  by  i|i^  =  >K£)j  that  is  we  have 


**  =  - 

G<|>  da 
,  n 

b 

(4.3) 


This  potential  tp*  is  known  in  principle  since  $  is  given  on  (b) .  A  choice 
of  methods  is  available  for  solving  the  nonhomogeneous  Fredholm  integral 
equation  (4.2).  In  particular,  a  natural  method  of  solution  consists  in  us¬ 
ing  an  iterative  procedure  based  on  an  appropriate  recurrence  relation.  An 
obvious  recurrence  relation  is  that  obtained  by  simply  replacing  <j>  by  <p  and 
^(k+1)  on  the  right  an(j  ieft  sides,  respectively,  of  equation  (4.2),  i.e. 


(l-w.)< 


(k+1) 


=  n  + 


<♦ 


(k)  A(k) 


)G  da 
n 


k>0 


(4.4) 


where  the  initial  (zeroth)  approximation  4/^must  be  specified  somehow.  The 

practical  usefulness  of  such  an  iterative  method  of  solution  depends  on  fast 

(k) 

convergence  of  the  successive  iterative  approximations  cp  ,  which  in  turn  is 
crucially  dependent  on  selection  of  a  "reasonably-good1'  initial  approximation 
In  some  particular  problems,  such  an  initial  approximation  is  readily 
available;  for  instance,  in  the  problems  of  wave  radiation  due  to  surging 
motions  of  a  slender  ship  or  to  heaving  and  pitching  motions  of  a  thin  (beam 
<<  length  and  draft)  ship,  we  may  simply  take  <ji^  =  0  (for  swaying  motions 
of  a  thin  ship,  on  the  other  hand,  this  simple  approximation  clearly  could 

not  be  used).  For  sufficiently  low  (or  high)  values  of  the  frequency  para- 

2  0  °° 
meter  f  =  w  L/g,  the  zero-  (or  infinite-)  frequency  potential  <p  (or  $  )  might 

also  be  used  as  an  initial  approximation;  more  generally,  a  relatively-simple 

interpolation  approximation  such  as  (<f>^+f(f>  )/(l+f),  which  clearly  reduces  to 

0  00 

<t>  and  41  as  f-+0  and  f-+<=°,  respectively,  would  seem  likely  to  provide  a 
reasonable  initial  approximation  4>^. 
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Equation  (4.2),  like  equation  (3.9)  of  course,  holds  not  only  for  £ 
on  the  body  surface  (b)+(c) ,  but  also  for  £  outside  (and  inside)  the  body. 

This  means  that  the  integral  equation  (4.2)  can  in  principle  be  used  to  determine 
the  potential  in  the  entire  solution  domain  (d)+(o)+(b)+(c) ,  specifically  by 
means  of  the  recurrence  relation  (4.4).  In  practice  however,  it  would  usually 
be  much  simpler  to  solve  for  <J>  on  the  body  surface  (b)  +  (c)  ,  and  —  in  the  event 
(rare  in  reality)  that  knowledge  of  <p  in  the  flow  domain  is  in  fact  required  — 
to  determine  <f>  outside  (b)  +  (c)  by  means  of  equation  (3.6a),  which  here  takes 
the  simplified  form 


G<(>  da  + 

n 


4>G  da 
n 


(4.5) 


Equation  (4.2)  also  holds  for  £  inside  the  body,  as  was  noted  previously. 

It  thus  might  appear  that  this  integral  equation  must  also  define  the  potential 
ip  inside  the  body.  This  result,  were  it  to  be  true,  would  certainly  be  quite 
surprising  (indeed  fundamentally  unacceptable)  for  it  would  mean  that  the  "exterior 
boundary-value  problem"  stated  by  equations  (1.6a,b,c,d)  would  define  a  solution 
in  the  "interior  domain"  (d^)+(o^)-(b)-(c) .  It  can  easily  be  shown,  however, 
that  equation  (4.2)  allows  the  potential  4>  to  be  extended  inside  (b)+(c)  in  an 
entirely  arbitrary  manner.  Indeed,  equation  (4.2)  can  be  written  in  the  form 


CP, 


G<f>  da  + 
n 


<j>G  da 
n 


(4.6) 


where  C  is  given  by  C  =  1-w.  +  /  G  da.  We  have  C  =  0  for  £  inside  (b)+(c) , 
as  may  easily  be  verified  (and  has  indeed  been  shown  in  the  previous  section) , 
so  that  equation  (4.6)  clearly  does  not  define  $  inside  (b)+(c) .  For  £  out¬ 
side  (b)+(c) ,  on  the  other  hand,  we  have  C  =  1,  and  equation  (4.6)  becomes 
equation  (4.5),  which  obviously  defines  ij>^. 

A  simple  explanation  for  the  origin  of  the  "waterplane  integral"  w^  = 
w(£;f,e,o^)  defined  by  formula  (3.9a)  in  the  integral  equation  (4.2)  in  the 
case  of  a  free-surface  piercing  body  [this  "waterplane  integral"  does  not 
appear  in  the  case  of  a  fully-submerged  body]  may  be  provided  by  considering 
a  free-surface  piercing  body  as  the  "zero-submergence  limit"  of  a  slightly 
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submerged  body  consisting  of  the  mean  wetted  body  surface  (b)  closed  bv  a 
horizontal  "lid",  (£)  say,  submerged  a  depth  6  below  the  mean  free-surface 
plane  z=0,  so  that  the  "lid"  (l)  becomes  the  "interior  waterplane"  GO  in 
the  limit  6=0.  For  the  fully-submerged  body  (b)+((!)  ,  the  integral  equation 
(4.2)  takes  the  form 


$*  =  -  G<f>  da  +  (<f>-4>+)G  da 

W  n  W  n 


=  -  |  G<)>  da  + 
>b  n 


f 

( 4>_ 4>*)  da  +  |  GfjOxdy  -  |  ( <j>-<(>*)  Gzdxdy  , 


where  the  fact  that  <j>  =  -ip  on  the  "lid"  (£)  was  used.  By  Expressing  the 

n  z  2  2  2 

terms  4>  and  G  in  the  forms  4>  -f(l+ie)  <t>  +  f(l+ie)  <J>  and  G  -f(l+ie)  G  + 
z  z  z  z 

o 

f(l+ic)  G,  we  may  then  obtain 


G<J>  da  + 
n 


(<M>*)G  da  +  4>.f  (1+ie)  2 
jb  *  n 


I  Gdxdy  + 

>1 


+ 


G [ 4>  -f  ( 1+ie)  ^4> ]dxdy  - 

>1  z 


( 4>— <J 


>L 


*) [Gz-f(l+ie)  G]dxdy 


2 

We  have  <f>  -f(H-ie)  <j>  =  0  on  the  mean  free-surface  plane  z=0,  so  that 
Z  2 

[<t>  -f(H-ie)  d>]  .  vanish  as  6->0.  Likewise,  we  have  ($-40  [G  -f  (1+ie)  ~G]->-0 

Z  Z=”0  *  z 

as  6-K).  The  last  two  integrals  on  the  right  side  of  the  above  equation  thus 
vanish  in  the  limit  6-H),  and  this  equation  may  then  be  seen  to  become  identi¬ 
cal  to  equation  (4.2)  in  this  limit.  Moreover,  the  term  w^<j>^  in  equation 
(4.2)  for  a  free-surface  piercing  body  (b)  can  be  seen  to  stem  from  the  effect 
of  the  "lid"  (£)  in  the  integral  equation  for  the  slightly  submerged  body 
(b )+(£). 

A  modified  form  of  the  integral  equation  (4.2)  can  be  obtained  by  expres¬ 
sing  the  potential  ^  =  <(>(£)  in  the  form  ^  =  k^x*,  where  x*  =  X(t)  is  some 
given  potential  (to  be  selected  somehow) ,  and  k^  =  k(£)  is  the  function  de¬ 
fined  by  the  relation  k^  =  <4*/x^«  By  substituting  this  expression  for  4>*  into 
equation  (4.2)  and  expressing  the  term  4>-<f’*  =  ky-k^x*  in  the  form  kvt(x~X*)  + 
(k-k^) x » we  may  obtain  the  following  integral  equation  for  the  function  k(£)t 


V 
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ru-w.)x*- 


(x-X*)Gda]kA  =  **  + 


(k--k*)xG  da 
*  n 


(A. 7a) 


The  last  integral  on  the  right  side  of  the  above  integral  equation  would 
evidently  vanish  if  k  k^.  for  x  and  r.  on  the  body  surface  (h)  ,  that  is  if 
k  =  <j>/x  is  constant  on  (b)  ,  or  if  x  Is  proportional  to  <p  on  (b) .  The 
integral  equation  (4.7a)  would  then  yield  the  explicit  solution 

k*  =  -  Jb(x-X*)Cnda]  ,  (4.7b) 


from  which  we  may  immediately  obtain  the  following  expression  for  the  velocity 
potential  4>*  =  k^x*: 


(X~X*)G  da] 
b 


(4.7c) 


An  alternative  form  of  equation  (4.7a)  may  be  obtained  by  multiplying  both 
sides  of  this  equation  by  x*  and  using  the  relation  4)  =  ky;  this  yields  the 
following  integral  equation  for  the  velocity  potential  <J> : 


ru-w*)x* 


(x-X*)GndaH*  =  x*<Jj*  +  (4>X*-<f>*x)Gnda 

b  b 


(4. 7d) 


which  thus  is  a  modified  form  of  the  original  integral  equation  (4.2).  The 
potential  defined  by  formula  (4.7c),  which  clearly  is  the  exact  solution  of 
the  modified  integral  equation  (4.7d)  if  x  is  proportional  to  over  the  body 
surface  (b) ,  would  then  seem  likely  to  provide  a  fairly-good  approximation  to 
the  exact  potential  if  x  is  so  chosen  that  k  =  $/x  does  not  vary  significantly 
over  (b).  An  obvious  choice  for  the  function  x  is  the  potential  defined  by 
formula  (4.3).  More  generally,  we  may  choose  x  as 


(4 . 7e) 

where  >.  is  some  (presumably  slowly-varying)  function  which  may  be  chosen  so 
as  to  render  the  potential  (4.7c)  as  "good"  an  approximation  as  possible  [it 
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may  readily  be  seen  from  equation  (4.7d)  that  the  potential  defined  by  formu¬ 
las  (4.7c)  and  (4.7e)  actually  becomes  exact  if  the  function  X  is  proportional 
to  <f>/i|/];  in  the  absence  of  any  information,  one  may  simply  take  X  =  1. 

By  substituting  expression  (4.7e)  for  x  into  expression  (4.7c)  for 
we  thus  may  obtain 

(4.8) 

Formula  (4.8),  where  the  potential  iJj  is  defined  by  formula  (4.3)  and  the 
free  auxiliary  function  X  may  simply  be  taken  as  X  =  1,  provides  an  explicit 
approximation  to  the  velocity  potential.  This  explicit  approximation  is  the 
simplest  result  obtained  in  the  above-described  theory  based  on  formulating 
an  integral  equation  for  the  velocity  potential. 
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5 .  The  Green  function:  single-integral  representations 

2  2  1/2 

It  is  convenient  to  introduce  the  notation  h  H  fp  =  f(x"  +y'  ) 

2  2  1/2 

v  =  fz',  and  d  =  (h  +v  )  =  fr' ,  so  that  we  have 

h  =  f  [  (x-C)2+(y-n)  2]1/2  s  a)2[(X-Xs)2+(Y-Ys)2]1/2/g  , 

v  =  f(z+?)  =  0)2(  Z+Zg)/g, 

d  i  f[(x-C)2+(y-n)2+(z+C)2]1/2  =  u2[(X-X  )2+(Y-Y  )2+(Z+Z  )2]1/2/g  , 

S  o  s 


where  (X,Y,Z)  and  (Xg,Yg,Zg)  are  the  dimensional  coordinates  of  the  field 
point  and  of  the  singularity,  respectively,  in  the  Green  function,  and  equa¬ 
tion  (1.3a)  was  used.  It  may  thus  be  seen  that  d  represents  the  adimensional 

2 

distance,  in  terms  of  L=g/u  as  reference  length,  between  the  field  point  and 
the  mirror  image  of  the  singularity  with  respect  to  the  free-surface  plane 
z=0 ,  while  h  is  the  horizontal  distance  (similarly  adimensional)  between 
these  two  points,  and  v  is  the  negative  of  the  vertical  distance. 

We  restrict  our  attention  to  the  limiting  case  e  =  +0  corresponding  to 
purely  oscillatory  flow.  By  performing  the  change  of  variable  k  =  fX  in  the 
inner  integral  in  expression  (2.10b),  we  can  express  the  Green  function 
G(x;^,f)  =  G(x; C, f ,e=+0)  in  the  form 


4irG(x;f,f)/f  =  -  1/fr  +  g(h,v)  , 


where  the  function  g(h,v)  is  defined  by  the  double  integral 


(5.1) 


g(h,v) 


1 

d 


2 

7T 


•tt  (•*>  (v-ihcos0)A 

d6  I  dA  e _ 

h  }0  A-(l+i0) 


(5.1a) 


By  performing  the  changes  of  variables  a  -  fp  and  B  =  fv  in  expression  (2.10a) 
we  may  obtain  the  alternative  double-integral  representation 


g(h,v) 


1 

d 


r 


,  2 ,  2,1/2  ,, 

v(p  +v  )  -ihp 

<„2-H,2)1/2-U+10> 


(5.1b) 


5  . 1  Haskind '  s_  in tc^r ja.  1_  represent a t i^o n 

We  start  by  expressing  the  double  integral  (5.1a)  in  the  form 

g(h ,  v)  =  -  1/d  -  (2/ti)  I  I(0;h,v)d0  ,  (  ^ .  2 

J0 

where  I(9;h,v)  is  the  inner  integral  given  by 


I  ( 6 ;  h ,  v) 


(v-ihcosO) X 

r 

e 

n  T-(l+iO) 


By  considering  the  contours  of  integration  in  the  complex  plane  X  =  X  +  i 
shown  in  figure  1  —  where  the  lower  and  upper  contours  are  selected  for 
0<9<tt/2  and  tt/2^0<tt,  respectively  —  we  can  express  the  integral  I  in  the 
forms 

-(hcosG+iv) t 

r 

I  =  I  - — - — : -  dt  for  0 <0 <tt /2  , 

J0 


I 


I 

1  0 


(hcos9+iv) t 
e _ 

+i 


v-ihcos0 

dt  +  2nie 


for  tt /2 <  0<tt 


Use  of  the  above  expressions  for  I  into  equation  (5.2)  then  readily 


yields 


tt/2  ~  (hcosO+iv)  t  (hcos9+iv)t 

i  2  ru  f  ,  ,.  2  (V  ( 


g(h, v)  =  -  -t  -  -  d9 

d  tt  ^  Jn  t-i 


0  J0 


dt - 

IT 


de 


it/ 2 


0 


t+i 


dt 


vz-tt  -ihcos9 
-4ie  e  d0  . 

tt  /  2 


After  performing  the  change  of  variable  0  =  in  the  last  two  integrals 
we  may  regroup  the  first  two  integrals  and  express  the  function  g(h,v)  in 
the  form 
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g(h,v)  =  W(h, v)  +  N(h,v), 


(5.3) 


where  the  functions  W(h,v)  and  N(h,v)  are  defined  as 


W(h,v)  =  -4ie 


tt/2  ihcosG 

e  d0  , 


tt/2 

N(h,v)  =  -l/d-(2/ir)  I  J(0;h,v)d0  , 

'o 


(5.4) 


with  J(0;h,v)  given  by 
-(hcos0+iv) t 


J  = 


r 


0 


r 


-(hcos0-iv) t 


t-i 


dt  + 


t+i 


dt 


(5.4a) 


The  integral  W(h,v)  may  be  expressed  in  terms  of  "standard  functions",  as 
may  be  seen  for  instance  from  equations  (9.1.18)  p.360  and  (12.1.7)  p.  496 
in  Abramowitz  and  Stegun  [8].  Specifically,  we  have 


W(h,v)  =  2Trexp(v)  [HQ(h)-iJ0  (h)  ] 


(5.5) 


where  and  are  the  usual  Struve  and  Bessel  functions,  respectively. 

By  performing  the  changes  of  variables  T=> (hcos0+iv) t  and  T=(hcos0-iv) t 
in  the  first  and  second  integrals,  respectively,  in  expression  (5.4a)  for  J, 
we  may  obtain 


r°° 


j  = 


-T 

e  dr 


-T 

e  dr 


Jq  T+v-ihcos0  Jq  T+v+ihcos0 


(5.6) 


By  grouping  these  two  integrals,  we  can  obtain 


J  =  2 


-T 

00 

e  ( i+v) di 

2  2  2 

0  (t+v)  +h  cos  0 
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Use  of  this  expression  for  J  into  equation  (5.4)  then  yields 


N(h,v) 


1 

d 


r<*  — T 

dTe  (x+v) 

0 


n/2 


d0 


2  2  2 
(i+v)  +h  cos  0 


where  an  interchange  in  the  order  of  integration  with  respect  to  t  and  0  was 
performed.  The  inner  Integral  (with  respect  to  0)  can  he  evaluated  (in  the 
usual  manner  by  transforming  this  integral  into  an  integral  around  the  unit 
circle  |z|  =  1  in  the  complex  plane  z  =  e1**)-  We  may  then  obtain 


N(h,v)  -  -  ±-  2 


-T 

on 

e  sgn(T+y)di  1 

0  [(r+v)2+h2]1/2  d 


2e 


“t 

oo 

e  sgnt  dt 

.2,.  2,1/2 

v  (t  +h  ) 


where  the  change  of  variable  t  =  t+v  was  performed.  The  last  expression  may 
readily  be  written  in  the  following  alternative  forms: 


N(h,v)  =  -  -  -  2e 


-t 


(t2+h2r1/2 


dt  +  2e 


(t2+h2)_1/2  dt 


N(h,v)  = 


-  -  -  4e 


-t 


(t2+h2)  1/2  dt 


+  2e 


-t 


2  2  -1/2 
(t  +h  )  '  dt 


The  first  integrals  in  the  above  alternative  expressions  can  be  expressed 
in  terms  of  "standard  functions",  as  may  be  seen  for  instance  from  equation 
(12.1.8)  p.  496  in  [8].  We  may  then  obtain 


N(h,v)  =  -  ~  +  ire  [Yq (h)-l^0(h)  ]  +  2 
.  v 

N(h,v)  =  -  i  +  2ire  [YQ (h) -SfQ (h)  ]  +  2 


"v  t+v  2  2  -1/2 
e  (tN-li  )  dt 

0 


"V  11 +V  2  2-1/2 
e  (t  +h  )  '  dt 


(5.7a) 


(5.7b) 


where  and  are  the  usual  Bessel  and  Struve  functions,  respectively,  as 
they  are  defined  in  [  8  ]  for  instance. 


A 
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By  using  the  above  alternative  expressions  for  N(h,v)  and  expression  .5.5) 
for  W(h,v)  in  equation  (5.3),  we  may  finally  obtain 


g (h , v)  =  Tre  [Y0(h)+H0(h)-2iJ0(h)]  -  i  +  2 


r-v  t+v 


,.2.,  2. -1/2 

(t  +h  )  dt. 


(5.8a) 


g(h,v)  =  2ne  [Y0(h)-iJg(h) ]  -  -  +  2 


t+v 


<tWr1/2  dt 


(5.8b) 


Expression  (5.8b)  is  identical  to  the  expression  obtained  by  Haskind  [1]  and 
given  in  Wehausen  and  Laitone  [9]  p.  477  eq.  (13.17').  The  modified  Haskind 
expression  (5.8a)  was  used  by  Yeung  [10]  p.  61  eq.  (4.5b),  and  was  also  recently 
rederived  by  Hearn  [11].  For  purposes  of  numerical  evaluation,  a  convenient 
alternative  form  of  the  integral  in  expression  (5.8a)  is  obtained  by  perform¬ 
ing  the  change  of  variable  x  =  -(t+v)/d.  This  yields 


I  g(h,v)  =  ire  [Yg(h)+HQ  (h)-2ijg(h)  ]  -  j  +  2 


a  -dr 


(l-2ax+T2)_1/2dT 


(5.8c) 


where  we  have  a  =  -v/d  by  definition  and  0<a<l.  This  modified  Haskind  integral 
representation  is  very  well  suited  for  evaluating  the  Green  function  for 
small  values  of  a.  Indeed,  for  v=0,  the  integral  in  formula  (5.8c)  vanishes, 
and  we  have  the  particularly  simple  expression 


g(h,v=0)  =  n[Y0(h)+ft0(h)-2iJ0(h)]-l/h 


(5.9) 


However,  Haskind's  integral  representation  is  clearly  not  well  suited 
for  evaluating  the  Green  function  for  values  of  a  close  to  1.  As  a  matter 
of  fact,  expressions  (5.8)  are  not  defined  for  a=l,  i.e.  for  h=0  and  v<0, 
so  that  these  expressions  can  only  be  used  for  v<0  and  h>0.  A  complementary 
single-integral  representation  that  is  well  suited  for  evaluating  the  Green 
function  for  v<0  and  small  values  of  h>0  will  now  be  derived. 
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5 . 2  The^  nearHF leld_in tegiral  r£p £es^ent at  i^on 

By  performing  the  changes  of  variables  t=T+v-ihcos0  and  t=T+v+ihcos0 
in  the  first  and  second  integrals,  respectively,  in  equation  (5.6),  we  may  obtain 

J  =  exp (v-ihcosO)Ej (v-ihcosO)+exp(v+ihcosO)Ej (v+iheosO)  » 

where  E1  is  the  usual  exponential  integral  defined  here  as  in  equation  (5.1.1) 
p.228  in  [8].  By  using  the  symmetry  relation  E^  (Z)  =  Ej  (Z) ,  given  by  equation 
(5.1.13)  p.229  in  [8]  for  instance,  we  may  then  obtain  J  =  2Re  exp (v+ihcosP) 
(v+ihcos0) .  Use  of  this  expression  for  J  into  equation  (5.4)  finally  yields 


N(h,v) 


-1/d  —  (A/tt) 


tt/2  Z 

Re  e  E.  (Z)d0 

0 


(5.10) 


where  Z  is  the  complex  number  defined  as  Z  =  v+ihcos0. 

By  using  formulas  (5.5)  and  (5.10)  in  formula  (5.3),  we  then  have 


'  g(h,v)  =  2 ire  [Hq (h) -iJQ (h)  ]  -  ^  ^ 


tt/2 


Re  e  E1(Z)d0;  Z  =  v+ihcos0  !  .  (5.11) 


For  the  sake  of  easy  reference  (and  for  reasons  which  will  become  clear  further 
on),  the  expression  for  the  Green  function  defined  by  formulas  (5.1)  and 
(5.11)  is  referred  to  as  the  "near-field  integral  representation"  of  the  Green 
function.  This  expression  has  also  been  obtained,  independently  and  in  a  dif¬ 
ferent  manner,  by  Martin  [2].  The  near-field  integral  representation  takes 
a  particularly  simple  form  for  h=0.  Indeed,  expression  (5.11)  then  becomes 


|  g(h=0, v)  =  1  / v  -2exp(v)  [ReE1(v+iO)+iir] 


(5.12) 


The  main  interest  of  the  integral  representation  (5.10),  in  comparison 
to  the  alternative  integral  representation  (5.7),  resides  in  that  this  expression 
can  be  used  to  obtain  an  ascending  series  useful  in  the  neighborhood  of  the  origin 
h=0=v,  i.e.  for  small  values  of  h  and  -v.  This  ascending  series  will  be  given 

in  section  7. 


40 


5.3.  The_  f^ar-_fiel(i  integral  rep  re  sen  tat.  i_on 

We  now  start  from  the  double-integral  representation  (5.1b),  which  we 
write  in  the  form 


g(h,v)  =  1/d  -  (2/tt ) 


oo 

0 


I (v;h,v) dv 


where  I(v;h,v)  is  the  inner  integral  defined  as 


I (v;h,v) 


’OO 


— oo 


,  2  2,1/2 

v(p  +u  )  -lhp 

e _ 

2  2  1/2 

(p  +v  )A  -(1+iO) 


(5.13) 


By  multiplying  the  numerator  and  denominator  of  the  integrand  of  the 

2  2  1/2 

inner  integral  I  by  the  expression  (p  +v  )  +  (1+iO) ,  and  by  rearranging 

the  denominator,  we  may  express  this  integral  in  the  form 


r 


r 

—  •» 


v(p2+v2)1/2-ihp 
[: +(l-v2+i0)1/2 


f(p2+y2)1/2  +  (1+iO) ] 
][p-(l-v2+10)1/2] 


dp  . 


2  1/2 

The  poles  ±(l-v  +i0)  of  the  integrand  of  the  above  integral  are  given  by 
21/2  21/2 

+  [(l-v  )  +i0]  if  0<v<l  and  +[0+i(v  -1)  ]  if  l<v<°°.  By  considering  the 

contour  of  integration  in  the  complex  plane  p  =  p  +  ip .  shown  in  figure  2, 

2  2  1/2  _  22  l/2r  1 

and  noting  that  we  have  (p  +v  )  =  -f-i(p^-v  )  for  p  =  +0+ip^  on  the 

two  sides  of  the  cut  defined  by  p^_  =  0  and  -  °°<p^<-v,  we  can  express  the 
integral  I  in  the  form 


.  .  .  ,  2  2,1/2 

hp ,+iv(p ,-v  ) 


,  .  ,  2  2,1/2 

-oo  hp  -iv(p  . -v  ) 


f-v  “h  ‘"h  v  '  ,-oo  npi-ivvpi-v  ) 

1  *  L  I(uj-V2)1/2  -  1  ld''1  +  L  .! 


id)j ^  -  2iriR. 


(5.14) 


In  this  expression  R  is  the  residue  at  the  pole 

2  1/2 

-i(p  -1)  if  l<v<°°,  so  that  we  have 


—  (1-v2)^2  if  0<v<l,  or  at 


41 


R(v;h ,v)  = 


■( 


-exp[v+ih(l-v2)1/2]/(l-v2)1/2  If  0<v<l 
iexp[v-h(  v2-l)  ]  /  (v^-l)  If  l<v<“> 


(5.15) 


2  •>]/■’ 

as  may  easily  be  found.  By  performing  the  changes  of  variables  t  =  (p.~\>~) 

091/2  1 

and  t=-(p.-v  )  in  the  first  and  second  integrals,  respectively, 

in  expression  (5.14),  we  can  finally  obtain 


I(v;h,v) 


.  ,  .2.  2.1/2  . 
r»  -h(t  +v  )  +ivt 


tdt  -  2iri  R(v;h,v)  , 


(5.16) 


where  R  is  given  by  equations  (5.15). 

By  using  equations  (5.16)  and  (5.15)  in  formula  (5.13),  we  may  then 
express  the  function  g(h,v)  in  the  form 


g(h,v)  =  W' (h,v)  +  N'(h,v)  , 


(5.17) 


where  the  functions  W'  and  N'  are  given  by 


W'  (h , v)  =  -4ie 


e _ 

/•i  2,1/2 

0  ( 1— v  ) 


dv  -4e 


-h(v2-l)1/2 

oo 

1  V-l)1/2 


dv 


(5.17a) 


N'  (h , v) 


1 

d 


2 

TT 


r 

dv 

J0 


r 

dp 


-h(p2+v2)1/2+ivp 
pe _ 

(p+i)(p2+v2)1/2 


(5.17b) 


The  two  integrals  in  equation  (5.17a)  can  be  expressed  in  terms  of  usual 

Bessel  functions,  as  may  be  seen  by  performing  the  changes  of  variables 

2  1/2 

v  =  sin9  and  v=(l+t  )  in  the  first  and  second  integrals,  respectively, 
and  by  using  equation  (9.1.18)  p.360  and  equations  (12.1.7)  and  (12.1.8) 
p.496  in  [8],  Specifically,  we  may  obtain 


W'(h,v) 


v 

2 ne  [YQ(h)-iJ0(h)] 


v 

-2irie  H 


(1) 

0 


(h) 


(5.18) 
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By  performing  the  changes  of  variables  |i=AcosO  and  v=AsinO  in  the 
double  integral  (5.17b),  we  may  express  the  function  N'(h,v)  in  the  form 


N'(h,v)  =  l/d-(2/u)  J(0;h,v)d0 

J0 


where  the  inner  integral  J  is  given  hy 


~(h-ivcos0) A 

r 

J(9;h,v)  =  -r - r -  AdA 

J0  A+isec0 


This  integral  may  be  expressed  as 


(5.1  ) 


’n-ivcos0 


-oo  -(h-ivcos0) A 

isecO  e _ 

■  0  A+isec0 


By  performing  the  change  of  variable  x=A+isec0,  we  may  then  obtain 


Z  /•“>  -(h-ivcos0)x  , 

T  isec9  .  dx 

J  =  — - —  -isec0e  e  — - 

isec0 


where  Z  is  the  complex  number  defined  as  Z  =  v+ihsec9.  The  change  of  variable 
t=(h-ivcos9)x  then  yields 


i  zr°°  Z 

J  =  isec9(—  -  e  dt)  =  isec0[l/Z-e  E1 (Z) ] 

Z  >Z 


Substitution  of  this  expression  for  the  inner  integral  J  into  formula  (5.19) 


then  gives 


tt  /  2  Z 


N'(h,v)  =  -1/d  -  (4/tt)  Im  e  E  (Z)sec0d0 

J0 


as  may  be  found  after  some  transformations. 


(5.20) 
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By  substituting  expressions  (5.18)  and  (5.20)  for  V  and  N',  respectively 
into  formula  (5.17),  we  can  finally  obtain 


;  V  1  A  "ILL* 

|  g(h,v)  =  2ire  [Y^hJ-iJ^h)  ]-  -  —  Ime  E^(Z)secBdO;  7.  =  v+ilisecO  | 


(5.21) 


This  expression  is  obviously  not  defined  for  h=0,  and  is  best  suited  for 
evaluating  the  function  g(h,v)  for  large  values  of  h.  The  expression  for  the 
Green  function  defined  by  formulas  (5.1)  and  (5.21)  will  thus  be  referred  to  as 
the  "far-field  integral  representation".  This  integral  representation  does  not 
seem  to  have  been  given  previously,  to  the  author's  knowledge.  There  is  a 
striking  similarity  in  form  between  the  far-field  representation  (5.21)  and  the 
near-field  representation  (5.11).  These  two  integral  representations  indeed 
are  complementary.  In  particular,  the  near-field  representation  (5.11) 
readily  provides  an  asymptotic  expansion  valid  for  large  values  of  d  and  small 
or  moderate  values  of  h/d,  that  is  in  a  sector  neighbouring  the  vertical 
axis  h=0,  while  an  asymptotic  expansion  valid  for  large  values  of  d  and  small 
or  moderate  values  of  -v/d,  that  is  in  a  sector  neighbouring  the  horizontal 
axis  v=0,  can  readily  be  obtained  from  the  far-field  representation  (5.21). 
These  complementary  asymptotic  expansions  are  given  in  the  following  section. 
Comparison  of  expressions  (5.8b)  and  (5.21)  show  that  these  expressions  are 
equivalent,  and  that  in  fact  we  must  have 


•tt/2  v+ihsec6  i-00  -dr  ~ 

Ime  E. (v+ihsec0)sec9d0  =  e  (1+2 jT+t  )~l/2  ^ 

Jo  J0  d 


as  may  be  obtained  by  performing  the  change  of  variable  r  =  -(t+v)/d  in  the 
integral  in  the  Haskind  expression  (5.8b) 
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6 .  The  Green  function:  asymptotic  expansions 

Let  us  first  consider  the  far-field  integral  representation  defined  by- 
expression  (5.21),  or  by  the  equivalent  equations  (5.17),  (5.18),  and  (5.20). 
We  define  the  integral  I^(h,v)  as 


(h,v) 


Im  (-2/tt) 


tt/2 

( l/Z) secOdf) ;  Z 

0 


v+ihsecO 


(6.1) 


We  have 


Ix(h,v)  =  1/d  ,  (6.1a) 

as  may  easily  be  verified,  and  indeed  was  already  used  in  the  derivation  of 
expression  (5.20),  By  using  equations  (6.1)  and  6.1a)  in  equation  (5.20),  we 
may  express  the  function  N' (h,v)  in  the  form 


N'  (h,v)  =  1/d  -(4/tt) 


(-tt/2  Z 

Im [e  E  (Z)  -1/Z]sec0d0;  Z 
0  1 


v+ihsec9 


(6.2) 


By  using  the  well-known  asymptotic  expansion 

n  n+1 

exp(Z)E  (Z)-l/Z  *v  £  (-1)  nl/Z  as  |  Z  |  -»-  °°,  with  lArgZl'-ir  ,  (6.3) 

n>l 

in  expression  (6.2),  we  may  obtain  the  asymptotic  expansion 


n 

N'(h,v)  1/d  +  2  £  (-1)  n(I  ,,(h,v)  a9  d-**,  with  h>0  ,  (6.4) 

.  n+i 


where  In+^(h,v)  is  the  integral  defined  by 


(-tt/2  n+1 

In+^(h,v)  =  Im(-2/7r)  j  (l/Z  )sec0d0;  Z  5  v+ihsec0 


It  may  be  seen  that  we  have  the  relation 
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In+l(h,V)  =  (_1/n)8In(h»v)/9v  « 


(6.5) 


from  which  we  may  obtain 


(-1)  n!ln+1(h,v)  =  3  I1(h,v)/3v  =  3  (l/d)/3v  ,  (n^l) 


(6 . 5a) 


where  equation  (6.1a)  was  used.  The  asymptotic  expansion  (6.4)  then  becomes 


N'  (h,v)  'v-  1/d  +  2  E  3  (l/d)/3v  as  d-*»,  with  h>0 
n>l 


(6.6) 


It  may  be  verified  that  we  have 


3  (l/d)/3v  =  P  (a)/d 

n 


(6.7) 


where  a  =  -v/d,  and  Pn(°0  is  a  polynomial  of  degree  n  in  a. 

By  using  equation  (5.18)  and  equations  (6.6)  and  (6.7)  in  equation  (5.17), 
we  may  finally  obtain 


g(h,v)  ^  2-rre  [Y.(h)-iJn(h) ]  +  1/d  +  2  1  P  (<x)/d  as  d^>  ' 
u  u  ,  n 

n>X 


(6.8) 


with  O^ar-v/d^l.  The  first  few  polynomials  Pn(«)  may  be  shown  to  be 


,  pl  =  a 


P2  =  -(l-3a  ) 


'  d  i2  n  5  2-» 

.  P  ^  =  — 3  a  (1-  -j  a  ) 


,2,.  ..  2  35  4. 
P^  =  3  (l-10a  +—  a  ) 


(6.8a) 


D  .3  2  14  2  21  4, 

P-  =  3  *5  a(l-~2-a  +-y<*  ) 


D  -2  ,.2,.  01  2  ,  .  4  231  6  . 

P,  =  -3  *5  (l-21a  +63a - ?-a  ) 

b  5  1 


We  now  consider  the  near-field  integral  representation  defined  by  express¬ 
ion  (5.11),  or  by  the  equivalent  formulas  (5.3),  (5.5),  and  (5.10).  We  define 
the  integral 
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■tt/2 

I  (h,v)  =  Re (—2/ tt)  (1/Z)d6  ;  2  ;  v+ihcos9  .  (6.9) 

0 

We  have 

I^h.v)  =  1/d,  (<>.9a) 

as  may  be  verified.  By  using  equations  (6.9)  and  (6.9a)  in  equation  (5.10''. 
we  may  express  the  function  N(h.v)  in  the  form 

fn  /2  z 

N(h,v)  =  1/d  -  (4/tt)  |  Re  [e  E  (Z)-l /Z  ]d0  ;  Z  =  v+ihcos0  .  (6.10) 

'o 

By  using  the  asymptotic  expansion  (6.3)  in  expression  (6.10),  we  may 
obtain  the  asymptotic  expansion 

n 

N(h,v)  1/d  +  2  £  (-1)  n!In+^(h,v)  as  d»,  with  v<0  ,  (6.11) 

n>l 

where  In+^(h,v)  is  the  integral 

■tt/2  n+1 

I  , .  (h ,  v)  =  Re  (— 2/tt)  (1/Z  )d6  ;  Z  =  v+ihcosB 

n+1  h 

Equations  (6.5)  and  (6.5a)  may  readily  be  verified  to  hold,  so  that  the 
asymptotic  expansion  (6.11)  becomes 

n  n 

N(h,v)  'v  1/d  +  2  £  3  (l/d)/9v  as  d->a>,  with  v<0  .  (6.12) 

n>l 

The  functions  N'(h,v)  and  N(h,v)  defined  by  equations  (6.2)  and  (6.10)  thus 
happen  to  have  the  same  asymptotic  expansion  as  d-*»,  although  these  expansions 
are  not  valid  in  the  same  regions  of  the  (h,v)  plane.  Specifically,  the  asymp¬ 
totic  expansions  (6.6)  and  (6.12)  are  not  valid  in  the  neighborhoods  of  the 


vertical  axis  h=0  and  of  the  horizontal  axis  v=0,  respectively.  It  may  be 

convenient  to  express  the  polynomials  P  (a)  in  equation  (6.7)  as  functions  of 

2  1/2  ^ 

3  =  h/d  [so  that  we  have  a  (1-3  )  ].  The  polynomials  P  («)  can  then  be 

n 

expressed  in  the  form 


Pn(a)  =  n!  Qn(3)  ,  (6.13) 


where  the  functions  Q  (3)  verify  Q  (0)=1,  since  we  have  P  (l)=n!  as  may  be  verified 

n  n  n 

from  equation  (6.8a). 

By  using  equation  (5.5)  and  equations  (6.12),  (6.7),  and  (6.13)  in  equation 
(5.3),  we  may  finally  obtain 


v  ^  n+1  ' 

g(h,v)  -v  2 ne  [H0(h)-iJQ(h)  ]  +  1/d  +  2  T.  nlQn(B)/d  as  d^>|  ,  (6.14) 

n>l  i 


with  0<^6=h/d<l.  The  first  few  functions  Q^(8)  are  given  by 


9  1/2 


=  (1-3  ) 


3  2 

q2  =  i  -je 


2  1/2  52 

Q3  =  (l-6z)  (1-|b^) 


!  %  ■  l-5«2+f  6‘ 


(6.14a) 


1/2 


,5  -  <1-62)‘"(l-762+f  e‘)  ,  q6  .  62+i|i6‘-  ^  [ 


The  difference,  6g(h,v)  say,  between  the  values  of  the  function  g(h,v) 
given  by  the  asymptotic  expansions  (6.8)  and  (6.14)  is  given  by 


5g(h,v)  =  2 ire  [YQ(h) -HQ(h)  ]  s  2tre 


-(l-e2)1/2d 


[Y0(6d)-Sf0(Bd)  ] 


(6.15) 


The  function  6g(h,v)  thus  is  exponentially  small  as  d-x»,  provided  we  have  0<8<1. 
It  may  then  be  seen  that  the  "transition  discontinuity"  6g(d,B=Bt)  due  to  the 


use  of  the  asymptotic  expansions  (6.14)  and  (6.8)  for  0<8<6t  and  respec¬ 


tively,  is  exponentially  small,  and  thus  is  negligible — in  an  asymptotic  sense 
in  comparison  with  the  algebraic  terms  l/d11  in  the  asymptotic  expansions.  An 
optimum  transition  between  the  asymptotic  expansions  (6.8)  and  (6.14)  may  be 
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determined  from  the  obvious  requirement  that  the  transition  discontinuity 
6g(d,8t>  is  a  minimum.  This  optimum  transition  then  is  given  by  the  solution 
of  the  equation  3  [6g(d,3) ]/33  =  0.  By  differentiating  equation  (6.15),  we 
may  then  obtain  the  following  equation  for  the  "transition  curve"  vf(h) 

-vt  =  h[H0(h)-Y0(h)]/[H1(h)-Y1(h)-2/Tr]  ,  (6.16) 

where  equations  (9.1.28)  p.  361  and  (12.1.11)  p.  496  in  [ 8  ]  were  used.  Tn 
particular,  equation  (6.16)  gives 


2  2  4 

-v  ^  h  (1+2/h  -30/h  +. . . )  as  h-x» 


(6.16a) 


as  may  be  obtained  by  using  equations  (12.1.30)  and  (12.1.31)  p.  497  in  [3  ]. 

By  substituting  expression  (6.16a)  into  equation  (6.15)  we  may  then  obtain 

2 

the  following  expression  for  the  transition  discontinuity  Sg^  (h)  ^  -4exp(-h  ) /h 

as  h-**>.  The  discontinuity  may  be  regarded  as  negligible  in  practice  if  it  is 

sufficiently  small  in  comparison  with  the  main  algebraic  term,  i.e.  1/d,  in 

2 

the  asymptotic  expansions  (6.8)  and  (6.14).  We  thus  require  that  4dexp(-h  ) /h 

be  smaller  than  the  desired  relative  accuracy,  e  say,  which  might  be  taken  as 

2  1/2  2 

e  =  .01  for  practical  applications.  This  then  yields  4(l+h  )  exp(-h  )<.01 

2  2  1/2  2  1/2  2 

[since  we  have  d  =  (h  +v  )  n,  h(l+h  )  as  h-»»  on  the  transition  curve  -v'^h  ], 

from  which  we  may  obtain  h>2.6  and  d>7.2. 

For  sufficiently  large  values  of  d  (say  for  d  greater  than  about  7 

according  to  the  foregoing  analysis),  the  function  g(h,v)  can  then  be  evaluated 

by  means  of  the  asymptotic  expression 


g(h,v)  ,v  W(h,v)  +  N(h,v)  as  d-*® 


The  function  W(h,v)  in  expression  (6.17)  is  given  by 


2-rexp(v)  [Y0(h)-iJQ(h)  ] 

v  (h) <v<0 

W(h,v)  = 

for 

2-exp(v) [HQ(h)-iJ0(h)  1 

— oo<v<V^  (h) 

(6.17) 


(6.17a) 
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where  the  transition  curve  vt(h)  is  defined  by  equation  (6.16).  The  function 
N(h,v)  in  expression  (6.17)  can  be  evaluated  by  means  of  either  one  of  the 
equivalent  asymptotic  expansions 


n+1 

1  N(h,v)  'v,  1/d  +  2  Z  P  (a)/d  r  1/d  + 
t  -  n 

I  n>l 


n+1 


2  Z  n.’Q  (6)/d 
,  n 
n>l 


as  d-x» 


(6.17b) 


where  a  =  -v/d  and  6  =  h/d.  and  the  functions  P  (a)  and  Q  (6)  are  given  by 
- 1 —  - — >  n  n 

equations  (6.8a)  and  (6.14a). 

The  error  associated  with  the  use  of  the  asymptotic  expansion  (6.17b) 

is  of  the  order  of  the  term  following  the  last  term  in  the  truncated  series 

(i.e.  the  first  discarded  term  in  the  series),  as  is  well  known.  The  requirement 

that  the  function  N(h,v)  be  evaluated  with  a  relative  accuracy  e  (say  with  e  =  .01 

in  practice)  may  then  be  approximately  expressed  by  the  condition  2 | Pn (a) | /dn<e , 

which  yields  d> [2 | P^(a) | /c ] The  function  d^(a;e)  =  [2 | P^Ca) I /e ] may  be 

evaluated,  notably  in  the  particular  cases  a=0  and  a=l .  In  the  particular 

case  a=l(6=0),  that  is  along  the  vertical  axis  h=0,  we  have  P  (l)=n!,  so  that 

X  /n  ^ 

we  may  obtain  d  (a=l,e)=(2n!/e)  .  For  e=.01,  we  may  then  obtain  d^=200,  d2=20, 

d^=10.63,  d^=8.32,  d,-  =  7.52,  d^=7.24,  d2=7.21,  dg=7 . 30 , . . . . In  the  particular  case 

a=0(8=l)  ,  that  is  along  the  horizontal  axis  v=0,  we  have  P?^  ^(0)=0  and  |P2n(0)|= 

l2 . 32. 52 . . . (2n-l) 2,  as  may  be  seen  from  equations  (6.8a).  We  may  then  obtain 

d2n(a=0,s)=[2.12.32.52.. . (2n-l)2/e]1/2n.  For  e=.01,  this  yields  d2=14.14,  d4=6.51, 

dg=5.96,  dg=6.21....  The  above  results  suggest  that  if  a  relative  accuracy  e=.01 

is  desired  (as  ought  to  be  sufficient  for  most  practical  applications),  it  may 

not  be  advantageous  to  use  more  than  the  first  five  terms,  i.e.  l^n<4,  in  the 

asymptotic  expansion  (6.17b);  furthermore,  it  appears  that  this  5-terms  asymptotic 

expansion  could  be  used  for  d  greater  than  about  8  if  a=l(h=0)  and  for  d  greater 

than  about  6.5  if  a=0(v=0).  These  values  of  d  fortunately  happen  to  be  about 

the  same  as  the  value  d=7.2  found  previously  from  the  requirement  that  the  transition 

discontinuity  in  the  value  of  the  function  W(h,v)  be  negligible.  On  the  basis 

of  the  foregoing  analysis,  it  may  thus  be  recommended  that  the  5-terms  asymototic 

expansion  (6.17b)  be  used  in  a  region  which  may  approximately  (and  tentatively)  be 

2  2 

defined  by  .the  equation  h  /45+v  /65  :>  1;  a  more  precise  numerical  determination  of 
the  domain  of  practical  usefulness  of  the  asymptotic  expansion  (6.17b)  is  of 
course  possible. 
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7 •  The  Green  function:  ascending  series 

In  this  section,  an  ascending  series  for  the  function  g(h,v)  is  obtained 
from  the  near-field  representation  given  by  formula  (5.11),  or  by  the  equivalent 
equations  (5.3),  (5.5)  and  (5.10).  Let  the  integrand  exp(Z)E^(Z)  in  equation 
(5.10)  be  expressed  in  the  form 

Z  Z  Z 

e  EX(Z)  =  -e  (£nZ+y)+e  ^ (Z)+£nZ+y]  .  (7.D 

Furthermore,  let  the  complex  number  Z  =  v+ihcos0  in  the  term  £nZ  be  written  in 
the  form 

_  d-v ,  2v  ,  .  2h 

Z  '  ~(d^+1d^r  COS0)  *  (7 •  2) 

tor  reasons  that  will  become  clear  further  on.  Also,  let  the  parameter  a  be 
defined  as 

a  =  h/ (d-v)  .  (7.3) 

2 

It  may  be  verified  that  we  have  0<o <1 ,  and  2v/ (d— v)  =  o  —1,  so  that  equation 
(7.2)  becomes 

d  v  2 

Z  =  — (  a  -l+i2crcos0) 

Use  of  this  expression  for  Z  in  equation  (7.1)  then  yields 

Z  v+ihcos0  d  2  Z 

e  E^(Z)  =  -e  [In  — j-+  y  +&n(a  -l+i2ocos0)]  +  e  [E^ (Z)+£nZ+y ] 

By  using  this  expression  for  the  integrand  exp(Z)E1(Z)  in  equation  (5.10),  we 
may  express  the  function  N(h,v)  in  the  form 


N(h,v)  =  -  i  +  2e  (£n^p+y)J0(h)  + 


2e  I  -  2J 


(7.4) 


5] 


where  equation  (9.1.18)  p.  360  in  [8]  was  used,  and  I  and  J  are  the  integral 


defined  as 


(•tt/2  ihcos0  2 

I  =  (2/tt) Re  I  e  £n[-(l-o  )+i2ocos0]d0  , 

Jo 

'TT  /  2  Z 

J  =  (2/tt)  Re  e  [E.  (Z)+^nZ+y  ]  d0  ;  Z  2  v+ihcosO 

Jo 


Substitution  of  expressions  (5.5)  and  (7.4)  for  the  functions  W  and  N  into 
equation  (5.3)  finally  yields  the  expression 


1/d  +  2e  [irH„(h)  +  (£n  V-' +Y-iir)j  (h)+Il 


The  ascending  series  for  the  above-defined  integrals  I  and  J  are  given  below. 
The  integral  J  is  considered  first. 

We  have 


Z  00  m  °°  k+1  k 

e  [E  (Z)+£nZ+y]  =  [  Z  Z  /m! ] [  Z  (-1)  Z  /k.k!]  , 

m=0  k=l 


as  readily  follows  from  the  ascending  series  of  the  functions  exp(Z)  and  E^(Z) 
The  above  product  of  series  may  be  expressed  in  the  form 


_  _  k+1  n  n 

M  _ DJ _ ]  z_  =  J  (  l  i)Z_ 

n=l  k=l  k  (n“k)!k!  n!  n=l  k=l  k)  n ! 


where  equation  (0.155,4)  p.4  in  Gradshteyn  and  Ryzhik  [12]  was  used.  The  well- 
known  binomial  expansion  formula  yields 


n  n  n  k  k  n-k  k 

Z  :  (v+ihcos0)  =  £  (P)i  h  v  cos  0 

k=0 


We  may  then  obtain 


n  n  k  2k  n-2k  2k 

ReZ  =  n!  E  (-1)  h  v  cos  0/ (2k) ! (n-2k) !  , 

k=0 


where  n'  is  defined  as  n'  =  n/2  if  n  is  even  or  as  n'  =  (n-l)/2  if  n  is  odd. 
We  thus  have 

Z  ■«>  n  n'  k  2k  n“2k  2k 

Ree  [E  (Z)+£nZ+Y]  =  E  (  E  ±)  [  E  (-1)  (2^,7n-2kY'  cos  0 1  • 

1  n=l  m=l  k=0  zk;  . 


Substitution  of  this  series  for  the  integrand  in  the  integral  (7.6)  and  term 
by  term  integration  finally  yields 


,  .  2k  n-2k 

,  !  ,  !  lw  ,  ,,  1.3.5.. .  (2k-l)  h  v 

J  =  v  +  l  (  E  —  )(v  /n!+  E  (-1)  7—7 - —  ToTTF? - oTT 

~  .m  ,  ,  2.4.6. ..(2k)  (2k)!(n-2k) 

n=Z  m=i  k=l 


n  1  _  / n/ 2 
!  ’  n  "  { (n-1)/; 


Use  of  equations  (7.10)  in  equation  (7.9)  then  yields 


2n 


2n-l 


1  l'°  +  n=l(_1)  l(2n)'  I2n  +  ‘  (2n~l )  !  i2n-l ! 


(7.11) 


The  integral  1^  can  be  expressed  as  a  contour  integral  around  the  unit  circle 
|z|  =1  in  the  complex  plane  z  =  exp(iO).  We  thus  have 


12  I’  = 

n  2tt 


M--(  z+-)(z-io)](z+-)  — 
z  o  z  z 


(7.12) 


2  =1 


By  using  the  binomial  theorem,  and  after  some  transformations,  we  may 


obtain 


,  ,  2n  n-1  2n-2k-l 

7<Z+I>  %fo(")<Z  +^2E«>  +  <Pf 


2n.  1 


i  i  2n-l  n-1  2n-2k-2 

hz*k  =  Z  (2in_1)  (z  +  ^ _ ) 


k~0 


2n-2k' 


Use  of  equations  (7.13a,b)  in  equation  (7.12)  then  yields 


in 


n-1 


Z2n  l(n  )X0  +  k^Q(k  )(12n-2k-l+  I2n-2k+l) 


^Vl  =  A<2r1)(I2n-2k-2  + 


k=0 


where  Im  and  l“  are  the  integrals  defined  by 


m 


C  "2T  j  £nfT’(z+^)(z"ic,)]z  dz,  m2  0  > 


(7.13a) 


(7.13b) 


(7.14a) 


(7.14b) 


(7.15a) 


j 
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I 

m 


_1_ 
2  n 


£n[—  (z  +—  )  (z-io)  ]  — 
z  a  m 

M=i 


m  >  2 


It  may  be  verified  that  we  have 


(7.15b) 


Furthermore,  we  have 

I_  =  I+  » 
m  m-2 


(7.16) 


(7.17) 


as  may  easily  be  verified  by  performing  the  change  of  variable  z=l/?  in  the 
integral  (7.15b).  Use  of  equations  (7.16)  and  (7.17)  in  equations  (7.14a,b) 
then  yields 


2n-l 

i2  *2n 


n-1 

£ 

k=0 


(2n)I+ 

^  k '  2n-2k- 


(7.18a) 


2n-2 

12  *2n-l 


n-1 

S  ( 

k=0 


2n-l 

k 


)I 


2n-2k-2 


(7.18b) 


m 

We  have  0 <a<l ,  so  that  the  function  £n[io(z+i/o) ]z  is  holomorphic  in  the 
region  |z|<_l,  and  the  integral  (7.15a)  becomes 


m 


1  f  o  /Z-ioN  m, 
yr  <-n(— : — )z  dz 

271  J|z|-1 


(7.19) 


This  integral  can  easily  be  evaluated;  we  have 

m  nri-1 

Im  1  i  o  / (m+1)  .  (7.20) 

Use  of  expression  (7.20)  for  the  integral  I+  in  equations  (7.l8a,b) 

m 

then  gives 
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n  2n  n-1  k  2n-2k 

(-D  I'  /(2n)!  =  (-2/2  )  Z  (-1)  o  /(2n-2k)k! (2n-k) ! 

k=0 


(7.21a) 


n  2n-l  n-1  k  2n-l-2k 

i (—1 )  II  i  / (2n-l)  !  =  (-2/2  )  Z  (-1)  o  / (2n-l-2k)k!  (2n-l-k)  !  .  (7.21b) 

k=0 


By  substituting  expressions  (7.16)  and  (7.21a,b)  for  the  integrals  1^ , 
1^  ^  in  equation  (7.11),  we  may  finally  obtain 


I  =  -2 


n 

Z 

n=l  k=0 


n-2k 


( 


(-1) 


k ! (n-k) !  n-2k 


)(£) 


h 

d-v 


n/2-1 

n'  =  { 

(n-l)/2 


l2n 


and 


(7.22) 


Equations  (7.22),  (7.8),  and  (7.7)  -where  the  classical  ascending  series  for 
the  Bessel  and  Struve  functions  -^(h)  and  ftg(h),  namely 


Vh> 


00  "'ll211  2  IT  ^  °°  11  211+1  I-)1)  2 

Z  (-D  (t  )  /(n!)  H  (h)  =  Z  (-1)  h  /I  .3  .5  ...(2n+l)  ,  (7.23a,b) 
n=0  1  n =0 


may  evidently  be  used  -provide  an  ascending  series  useful  for  evaluating  the 
function  g(h,v)  for  small  and  moderate  values  of  d. 
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8 .  The  Green  function:  one-dimensional  Taylor  series  expansions 

The  near-field  representation  (5.11)  may  be  used  to  obtain  a  Taylor 
series  expansion  of  the  function  g(h,v)  in  the  neighborhood  of  the  vertical 
axis  h=0.  Let  I(h,v)  be  the  function  defined  as 

rv/2  Z 

I(h,v)  =  Re ( 4 / it )  e  E1  (Z)d0;  Z  v+ihcos0  .  (8.1) 

J0 

By  expanding  the  function  I(h,v)  in  a  Taylor  series  about  tiie  axis  h=0,  we  may 
obtain 


°°  n  n 


I (h , v)  =  E  h  [3  I(h,v)/3h  ],_n/n'.  . 
n=0  h'° 


(8.2) 


Differentiation  of  both  sides  of  equation  (8.1)  yields 


3  I(h,v)/3h  =  Re(4/7t) 


tt/2  n  Z 


n  n  n 


[d  e  E^(Z)/dZ  ]i  cos  0d0. 


(8.3) 


From  the  definition  of  the  exponential  integral  function  E^(Z),  given  for 
instance  by  equation  (5.1.1)  p.  228  in  [8],  it  may  be  verified  that  we  have 


n  Z 


n  Z 


d  e  E  (Z)/dZ  =  e  E  (Z)  +  Z  (k-l).'/(-Z)  for  n_>l. 
1  k=l 


(8.4) 


By  using  equation  (8.4)  in  equation  (8.3),  we  may  obtain 


rni  1 


3h 


2n 


■  2<-i)n  2:lt.::<2n)1) [*v<“*1(v+io> + fx- 1  *>■■  ^  <8-5a> 

i=0  k=l  (-v) 


while  for  n=0  we  have 


I (h=0,v)  =  2e  ReE1(v+iO) 


(8.5b) 
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We  have  ReiE^ (v+iO)  ImE^ (v+iO)  — tt,  as  may  be  found  from  the  ascending  series 
for  the  exponential  integral  given,  for  instance,  by  equation  (5.1.11)  p.  229 
in  [8],  It  may  then  be  seen  that  we  have 

v 

3I/3h|  =  4e  ,  (8.5c) 

h=0 


and 


a2n+1i 


3h 


2n+l 


1  h=0 


-  2.4,6... (2n)  v  . 

3.5.7. .. (2n+l)  6  for  n-1‘ 


(8. 5d) 


Use  of  equations  (8.5a,b,c,d)  in  the  series  (8.2)  then  yields 


2n 


Kh.v)  -  2e  EeEl(v+iOm+  E  (-1)  +  + 


n=l 


2n+l 
1 

n=1  .  (2n+l)  (2n+l)1 

where  the  function  J(h,v)  is  defined  by  the  series 


v  n 

+  4e  Th+  £  (-It  2.4.6..  •  (2n)  h _  . 

_ V  }  3.5.7...(2n+l)(2n+l)!  J 


(8.6) 


00  n  2n 
J(h,v)  =  £  (-1)  [  £ 


2n 


(k-l)K  1.3.5. ..(2n-l)  h 


n=l  k=l  (-v) 


k  2.4.6.. .(2n)  (2n) ! 


(8.6a) 


The  two  series  between  brackets  in  equation  (8.6)  are  the  ascending  series  for 
the  functions  JqOO  and  (tt/2)Hq(Ii),  respectively,  as  may  readily  be  verified 
from  equations  (7.23a,b),  so  that  equation  (8.6)  becomes 


Kh.v)  =  2e  [ReE1(v+iO)J0(h)+TrH0(h)]  +  2j(h,v)  . 


(8.7) 


The  series  (8.6a)  for  the  function  J(h,v)  may  be  expressed  in  the  form 
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J(h,v) 


00  n  2n-l 

£  (-D  [  £ 

n=l  m=0 


(2n-l-m) ! 


(-v) 


2n-m 


2n 

1.3.5..  . (2n-l)  _h _ 

2.4.6. . .(2n)  (2n) ! 


(8.7a) 


By  substituting  expression  (8.7)  for  the  integral  (8.1)  into  formula 
(5.11),  we  may  finally  express  the  function  g(h,v)  in  the  form 

(8.8) 

The  series  (8.7a)  for  the  function  J(h,v)  can  be  written  in  the  form 

(8.8a) 

where  P  (v)  is  the  polynomial  of  degree  2n-l  given  by 

(8.8b) 

rtTe  have  J(h=0,v)=  0  and  Jp(h=0)=l,  so  that  expression  (8.8)  becomes  identical 
to  expression  (5.12)  in  the  limit  h=0. 

A  complementary  Taylor  series  of  the  function  g(h,v)  in  the  neighborhood 
of  the  horizontal  axis  v=0  can  be  obtained  from  the  Haskind  integral  representation 
(5.8a).  Let  l(h,v)  be  the  function  defined  as 


g(h,v)  =  -l/d-2JQ(h)exp(v)  [ReE1(v+iO)+iTr]-2J(h,v) 


rv  t+v  2  2-1/2 

I(h,v)  =  I  e  (t  +h  )  dt  .  (8.9) 

J0 


Expansion  of  this  function  in  a  Taylor  series  about  the  axis  v=0  yields 


I(h,v) 


00  n  n 


n 


E  v  [3  I (h ,v) / 3v  ]v=Q/n! 


(8.10) 


where  the  fact  that  l(h,v=0)=0  was  used.  It  can  be  verified  that  we  have 
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2n+l 

a _ i 

.  2n+l 

9v 


v=0 


n  k+1 

z  (-D 

k=0 


2  2  2  2 

1  .3  .5  .. . (2k-l) 

,  2k+l 
h 


2n+2 

9 _ I 

2n+2 

9v 


for  n>0. 


v=0 


Use  of  equation  (8.11)  in  equation  (8.10) 


2n+l 


I (h, v)  =  E 


k+1 


n=0(2n+1>! 


-d+ 


2n+2 


)[  Z  (-1) 
k=0 


yields  the  series 


2  2  2  2 

1  .3  .5  . . . (2k-l ) 

,  2k+l 
h 


which  may  be  written  in  the  equivalent  form 


I (h,v) 


E 

n=0 


2n+l 

(-v) 
(2n+l)  ! 


(1+  2^>t  1  (-1) 


n-m 


m=0 


2  2  2  2 

1  .3  .5  ... (2n-l-2m)  . 

,  2n+l-2m 
h 


By  using  equation  (8.9)  in  equation  (5.8),  we  may  then  express 
function  g(h,v)  in  the  form 

g(h,v)  =  -1/d  +  7rexp(v)[Y0(h)+H0(h)-2iJ0(h)J  +  2I(b,v) 

The  series  (8.12)  for  the  function  I(h,v)  may  be  written  as 


n 


I(h,v)  =  ^(l+f)  +  I 
n=l 


1.3.5..  ,(2n-l)  .  -v  . 

2.4.6. . .(2n)  (  h  ; 


2n+l 


(1+ 


2S«>!’„<h) 


where  Pn(h)  is  the  polynomial  of  degree  2n  given  by 


P  (h) 
n 


1  + 


n  m  2m  2  2  2  2 

E  (-1)  h  / (2n-l)  (2n-3)  (2n-5)  . . . (2n+l-2m) 
m=l 


.11) 


(8.12) 

the 

(P.13) 

(8.13a) 

(8.13b) 


In  the  limit  v=0,  we  have  1=0,  and  expression  (8.13)  becomes  identical  to 
expression  (5.9). 


9 .  The  gradient  of  the  Green  function 

The  vertical  derivative,  Gz,  of  the  Green  function  G(x;£,f)  may  be 
expressed  directly  in  terms  of  the  function  g(h,v),  as  will  now  be  shown.  Following 
an  idea  used  by  Eggers  [4]  for  the  similar  problem  of  ship  wave  resistance,  we 
express  the  Green  function  in  the  alternative  forms 

4ttG(x;4»  f)  =  -l/r  +  l/r'+H+(p,z';f)  =  -1/ r-l/r  '+H~  (p  ,z '  ;  f )  ,  (9.1a,b) 

2  21/2  •)  2  1/’ 
where  we  have  p  =  [(x-£)  +(y-n)  1  ,z'  =  z+4  ,  r  =  [p  ■ +(z-r,)  ]  “  , 

2  2  1/2  2 

r'  =  (p  +z'  )  ,  and  f  =  m  L/g,  as  was  defined  previously.  Although  the  precise 

expressions  for  the  function  H+(p,z';f)  and  H  (p,z';f)  can  evidently  be  readily 
obtained  from  the  analysis  in  the  previous  sections,  for  instance  by  setting 
e  =  +  0  in  expressions  (2.10a)  and  (2.10b),  the  precise  forms  of  these  functions 
is  actually  not  required  here.  By  using  equations  (9.1a,b),  we  may  obtain 

4ir(G  -fG)  =  (-1/r-l/r')  +  H~  -  f (-1/r+l/r')  -  fH+  .  (9.2) 

The  free-surface  condition  (2.13b)  shows  that  we  have  G  -fG=0  on  z=0  if  C<0.  It 

z 

may  also  readily  be  seen  that  -1/r+l/r' =0  on  z=0  and  (1/r+l/r')  =0  on  z=0  if  ;<0. 

^  z 

It  then  follows  from  equation  (9.2)  that  we  have  H  -fH  =0  on  z=0  if  £<0.  This 

z  + 

relation  however  must  hold  for  z<0  since  the  functions  H  and  H  depend  on  z+i;. 

We  thus  have  H^-fH+sO,  as  may  also  readily  be  verified  From  equations  (2.9a) 
and  (2.9b).  Equation  (9.2)  then  becomes 


4t (G  -fG)  =  -(1/r+l/r')  +  f(l/r-l/r') 

z  z 

By  using  expression  (5.1)  for  the  Green  function,  that  is 


4’rG(x;  £ ,  f )  / f  =  -1/fr  +  g(h,v) 


in  equation  (9.3),  we  may  finally  obtain 


4TrC>z/f2  =  (z-;)/f2r3  +  v/d3  -  1/d  +  g(h,v)  , 


(9.3) 


(9.4) 


(9.5) 


where  v  :  fz',  h  =  fp,  and  d  =  fr' ,  as  was  defined  previously.  Expression  (9.5) 
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for  he  vertical  derivative  G z  of  the  Green  function  was  obtained  previously 
by  Martin  [2],  in  a  different  manner.  The  horizontal  derivatives  G  and  G,  of 
the  Green  function  G  may  readily  be  obtained  by  differentiating  expression  (9.4). 
Specifically,  we  may  obtain 


with  G„  =  3G/3p  given  be 

P 


(9. 6a, b) 
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